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Chapter 1 

Introduction and Overview 



^ 1% ^Mw-M^ #r^ ^5?^ wTFf ycj^ci | 



Charak Samhita (First book of Aurveda) 

Even if we separately know all the organs making up a human body, we still fall short of knowing the 

whole body. 



In this chapter we introduce the research areas dealt with in this thesis 
and give some background material. 

1.1 Quantum entanglement 

Entanglement is a subtle and eluding property of quantum systems com- 
prising many parts. Entanglement induces correlations between the mea- 
surable properties of different parts of a quantum system which cannot be 
reproduced by any procedure involving only the local operations (LO) and 
classical communication (CC) between various parts of the system jl]. In 
consonance with this, entanglement in a quantum system cannot increase 
(or be created) via LOCC. This principle is connected to another intriguing 
property of entanglement: a multipartite quantum system can get entan- 
gled in various inequivalent ways, which cannot be transformed into each 
other via LOCC. However, the most challenging aspect of entanglement is 
that it cannot be 'built in parts', that is, the entanglement of N parts is 
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not a sum or a simple function of the entanglement of M(< A^) partite 
subsystems f2]. 

The concept of entanglement has played a crucial role in the devel- 
opment of quantum physics. In the early days entanglement was mainly 
perceived as the qualitative feature of quantum theory that most strikingly 
distinguishes it from our classical intuition. The subsequent development of 
Bell inequalities made this distinction quantitative, and therefore rendered 
the nonlocal features of quantum theory accessible to experimental verifi- 
cation P, HI, [3] . Bell inequalities may indeed be viewed as an early attempt 
to quantify quantum correlations that are responsible for the counterin- 
tuitive features of quantum mechanically entangled states. At the time 
it was almost unimaginable that such quantum correlations between dis- 
tinct quantum systems could be created in well controlled environments. 
However, the technological progress of the last few decades means that we 
are now able to coherently prepare, manipulate, and measure individual 
quantum systems, as well as create controllable quantum correlations. In 
parallel with these developments, quantum correlations have come to be 
recognized as a novel resource that may be used to perform tasks that are 
either impossible or very inefficient in the classical realm. These devel- 
opments have provided the seed for the development of modern quantum 
information science. 

Given the status of entanglement as a resource it is quite natural and 
important to discover the mathematical structures underlying its theoreti- 
cal description. We will see that such a description aims to provide answers 
to three questions about entanglement, namely (1) its detection and clas- 
sification, (2) its creation and manipulation and, (3) its quantification. 

In this thesis, we deal with the second and the third problem. We 
have used the geometric (Bloch representation) approach for studying the 
creation and quantification of entanglement we give a geometric measure 
for quantifying the entanglement of multipartite pure states of fermionic 
systems. Our measure satisfies all properties for a good measure of entan- 
glement. 

We have loosely described entanglement as the quantum correlations 
that can occur in many-party quantum states. This leads to the question 
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what differentiates quantum correlations from as implied by entanglement 
classical correlations? In the context of quantum information one of the 
precise way to define classical correlations is via LOCC operations. Clas- 
sical correlations can be defined as those that can be generated by LOCC 
operations. If we observe a quantum system and find correlations that can- 
not be simulated classically, then we usually attribute them to quantum 
effects, and hence label them quantum correlations. The entanglement is a 
resource because it hfts the so-called LOCC constraint, i.e. entanglement 
and LOCC together can perform tasks that cannot be accomplished by 
LOCC alone. Using LOCC-operations as the only other tool, the inherent 
quantum correlations of entanglement are required to implement general, 
and therefore nonlocal, quantum operations on two or more parts P, E] . As 
LOCC-operations alone are insufficient to achieve these transformations, 
we conclude that entanglement may be defined as the sort of correlations 
that may not be created by LOCC alone. 

Entanglement has proved to be a vital physical resource for various kinds 
of quantum-information processing, including quantum state teleportation 
pi?l [TT] , cryptographic key distribution |T2] , classical communication over 
quantum channels |[I3l [HI [15], quantum error correction [il6j, quantum 
computational speedups IE] , and distributed computation ||18l [I9] . Fur- 
ther, entanglement is expected to play a crucial role in the many particle 
phenomena such as quantum phase transitions, transfer of information 
across a spin chain [TQl [21] etc. Therefore, quantification of entanglement 
of multipartite quantum states is fundamental to the whole field of quan- 
tum information and in general, to the physics of multicomponent quantum 
systems. 

Whereas the entanglement in pure bipartite state is well understood, 
the understanding of entanglement in mixed bipartite state is far from 
complete. In section 1.1, we review the entanglement of bipartite quantum 
system. We will state the available measures and criteria for detecting 
entanglement for both bipartite pure and mixed states. In section 1.2, 
we deal with multipartite entangled states. In section 1.3, we explain the 
capability of creating entanglement for a general physical interaction acting 
on two qubits. In section 1.4, we discuss entanglement in indistinguishable 
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particle systems. In section 1.5, we summarize the quantum discord as a 
measure of "quantumness" of the system. The material in section 1.3, 1.4, 
1.5 and 1.6 forms a background for chapters 2, 3, 4, and 5. Section 1.7 is 
chapterwise summary. 



1.2 Bipartite Entanglement 

In this section, we define the entanglement in bipartite quantum states. 
We review the work that has been done in the bipartite systems. Consider 
a system consisting of two subsystems. Quantum mechanics associates to 
each subsystem a Hilbert space. Let 1-La and I-Lb denote these two Hilbert 
spaces, let \i)A (where i = 1,2,3, •••) represent a complete orthonormal 
basis for Ha and \j)b (where j = 1,2,3, •••) a complete orthonormal 
basis for TIb. Quantum mechanics associates with the system, i.e. the 
two subsystems taken together, the Hilbert space spanned by the states 
\i)A \J)b- In the following, we will drop the tensor product symbol <S> 
and write \i)A'S> \j)b as |^)a|J)b, and so on. Any linear combination of the 
basis states \i)A\j)B is a state of the system, and any state |^)as of the 



system can be written [22 



ab = ^c^j\i)a\j)b, (1-1) 



^J 



where the Qj are complex coefficients, we take |^)as to be normahzed, 
hence J2ij hj? = '^■ 

If we can write \t/j)AB = \iP^^'')a\iI^^^^)b, we say the \t/j)AB is product 
state (separable state). If \iIj)ab is not a product state, we say that it is 
entangled. 

By using local operators and classical communication (LOCC) any state 
\i^)AB of two subsystems A and B can be transformed to the form [23, El 



AB = Y.d,\(l),)A\(l)[)B; k<dim{nA^nB), (1.2) 
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where the positive coefficients di are called Schimdt coefficients. The 
state is entangled if at least two coefficients do not vanish. Pure entangled 
state contains quantum correlation which can not be simulated by any 
classical tools. A fundamental Theorem was proved by Bell [S], who showed 
that if the constraint of locality was imposed on the hidden variables, then 
there was an upper bound on the correlations of results of measurements 
that could be performed on the two distant systems. That upper bound, 
mathematically expressed by Bell's inequality j3], is violated by some state 
in quantum mechanics, thus the state contains quantum correlation which 
is Non-local property of quantum state [21] . 

However, in real conditions, owing to interaction with the environment, 
called decoherence, we encounter mixed states rather than pure ones. A 
mixed state is a classical mixture of pure quantum states [22] • These 
mixed states can still possess some residual entanglement. A mixed state 
is considered to be entangled if it is not a mixture of product states [25] . 
In mixed states the quantum correlations are weakened and hence the 
manifestations of mixed state entanglement can be very subtle [2S] • From 
the definition of entanglement of mixed state it is difficult to apply this 
definition directly to know the quantum state is entangled or not, because 
the mixed state contains both classical and quantum correlations, and can 
be prepared using infinite possible ensembles. 

For pure states, it is easily shown that the CHSH inequality is violated 
by any nonfactorable state [27[ EB] , while on the other hand a factorable 
state trivially admits a (contextual) LHV model [4] . 

For mixed states, Werner [25] constructed a density matrix p^j for a pair 
of spin-j particles. Werner's state p^ can not be written as a sum of direct 
products of density matrices, ^jCjp^ pf, where A and B refer to the 
two distant particles and j runs over the states in the ensemble. There- 
fore, genuinely quantum correlations are involved in p^j. Nevertheless, for 
any pair of ideal local measurements performed on the two particles, the 
correlations derived from p^ not only satisfy the CHSH inequality, but, as 
Werner showed [25] , it is possible to introduce an explicit LHV model that 
correctly reproduces all the observable correlations for these ideal measure- 
ments [21]. Thus for mixed states entanglement and nonlocality are two 
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different resources. 



1.3 Multipartite Entanglement 

Multiparticle entanglement is genuinely different from entanglement in 
quantum systems consisting of two parts. To understand what is so dif- 
ferent consider, say, a quantum system that is composed of three qubits. 
Each of the qubits is to be held by one of three laboratories distantly 
separated. It may come as quite a surprise that states of such composite 
quantum systems may contain tripartite entanglement, while at the same 
time showing no bi-partite entanglement at all. In contrast to the bipartite 
setting, there is no longer a natural "unit" of entanglement, the role that 
was taken by the maximally entangled state of a system of two qubits. 
Quite strikingly, the very concept of being maximally entangled becomes 
void. Instead, we will see that in two ways there are "inequivalent kinds 
of entanglement" . Consider multipartite entanglement of pure quantum 
states. A theory of entanglement should not discriminate states that dif- 
fer only by a local operation. Here, "local operation" can mean merely 
a change of local bases (LU operations) or, else, general local quantum 
operations assisted by classical communication, that are either required to 
be successful at each instance (LOCC) or just stochastically (SLOCC). For 
each notion of locality, local unitary operation (LU) or Local operations 
and classical communication (LOCC) or just stochastic-LOCC (SLOCC), 
the questions that have to be addressed are how many equivalence classes 
exist, how are they parameterized and how can one decide whether two 
given states belong to the same class? 

For the case of bi-partite qubit states, two quantum states are LU- 
equivalent if and only if their respective Schmidt normal forms coincide. 
All classes are parameterized by only one real parameter. Some simple 
parameter counting arguments show that in the case of N-qubit systems 
the situation must be vastly more complex. Indeed, disregarding a global 
phase, it takes 2^"*^^ — 2 real parameters to fix a normalized quantum state 
in H = (C^)*^^. The group of local unitary transformations SU{2) x • • • x 



1.3 Multipartite Entanglement 



SU{2) on the other hand has 3N real parameters [2]. Therefore, one needs 
at least 2^+^ — 3N — 2 real numbers to parameterize the sets of inequivalent 
pure quantum states [29j. This lower bound turns out to be tight [30]. It 
is a striking result that the ratio of non-local to local parameters grows 
exponentially in the number of systems. In particular, the finding rules out 
all hopes of a generalization of the Schmidt normal form. A general pure 
tripartite qubit state, say, cannot be cast into the form sin9\()00)-\-cos9\lll) 
by the action of local unitaries ||5] • Considerable effort has been undertaken 
to describe the structure of LU- equivalence classes by the use of invariants 
or normal forms [291 EO, [3T], [321 E3]- Acm et al. [34j have proved for any 
pure three-qubit state the existence of local bases which allow one to build 
a set of five orthogonal product states in terms of which the state can be 
written in a unique form. This leads to a canonical form which generalizes 
the two-qubit Schmidt decomposition. It is uniquely characterized by the 
five entanglement parameters. When one deals with SLOCC operations, 
the group of SLOCC, SL^C^) x • • • x SLiC^) has 2^+i-6iV- 2 parameters 
that are necessary to label SLOCC equivalence classes of qubit systems. It 
turns out that the three-qubit pure states are partitioned into a total of Six 
SLOCC-equivalent classes [25] • The picture is complete for three-qubit : 
any fully entangled state is SLOCC-equivalent to either \GHZ) or \W) [[35]. 
Three-qubit W-states and GHZ-states have already been experimentally 
realized, both purely optically using postselection [Ml EZ] and in ion traps 
[SB]- The two states behave differently, however, if a system is traced 
out. Specifically, tracing out the first qubit of the GHZ state will leave the 
remaining systems in a complete mixture. For \W) will leave the remaining 
systems in a mixed entangled bipartite state. Thus, the entanglement 
of \W) is more robust under particle loss than the one of \GHZ) [33]. 
From point of view of asymptotic manipulation of multipartite quantum 
states, there is no longer a single essential ingredient as in bipartite the 
maximally entangled state or EPR-state, but many different ones. In the 
multi-particle case, however, it is meaningful to introduce the concept of a 
minimal reversible entanglement generating set (MREGS). An MREGS 
(S is a set of pure states such that any other state can be generated from 
S by means of reversible asymptotic LOCC. It must be minimal in the 
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sense that no set of smaller cardinality possesses the same property 
401 I4T] . Yet, it can be shown that merely to consider maximally entangled 
qubit pairs is not sufficient to construct an MREGS f40]. To find general 
means for constructing MREGS constitutes one of the challenging open 
problems of the field: as long as this question is generally unresolved, the 
development of a "theory of multi-particle entanglement" in the same way 
as in the bi-partite setting seems unfeasible. 



1.4 Generation of entangled quantum systems : en- 
tanglement capacity 

Since entanglement cannot be increased (generated) via LOCC, the in- 
teraction between various parts of a quantum system is needed to gen- 
erate and increase its entanglement. Thus a quantum system evolves 
to generate entanglement provided its parts interacts. Any Hamiltonian 
'Hab 7^ Ha+T-Lb that is not a sum of local terms couples the systems A and 
B. Together with local operations, the coupling can be used to generate 
entanglement [ H21 H51 IH] , to transmit classical and quantum information 
151 1151 1151 H7] , and more generally, to simulate the bipartite dynamics of 



some other Hamiltonian Hab and thus to perform arbitrary unitary gates 
on the composite space 1-Lab = 'Ha ® ^b [IHl US, SS] • Much experimental 
effort has been devoted to creating entangled states of quantum systems, 
including those in quantum optics, nuclear magnetic resonance, and con- 
densed matter physics [5D]. Determining the abifity of a system to create 
entangled states provides a benchmark of the "quantumness" of the sys- 
tem. Furthermore, such states could ultimately be put to practical use in 
various quantum information processing tasks, such as superdense coding 
fSIJor quantum teleportation |55] . 

The theory of optimal entanglement generation can be approached in 
different ways. For example, Ref. [12] considers single-shot capacities. In 
the case of two-qubit interactions, and assuming that ancillary systems are 
not available, Ref. [12] presents a closed form expression for the entan- 
glement capacity and optimal protocols by which it can be achieved. In 
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contrast, Ref. [13] considers the asymptotic entanglement capacity, allow- 
ing the use of ancillary systems, and shows that when ancillas are allowed, 
the single-shot and asymptotic capacities are in fact the same. However, 
such capacities can be difficult to calculate because the ancillary systems 
may be arbitrarily large. 

In the past few years there has been a considerable increase in exper- 
imental activity aiming to create entangled quantum states. One reason 
is the potential applications of entanglement to quantum information pro- 
cessing. Creating entanglement has been possible in quantum optics for 
more than a decade; however, now many new communities, working in a va- 
riety of experimental areas (for example, NMR, condensed matter physics) 
are also joining the field fW\. In general, entanglement between two sys- 
tems can be generated if they interact in a controlled way. However, in 
most experiments these interactions are weak which makes the production 
of entanglement a very difficult task. Thus, it would be very convenient 
to have a theory which would provide us with the best way of exploiting 
interactions to produce entanglement, we try to analyze the entanglement 
capabilities of Hamiltonians. In particular, we would like to answer ques- 
tions such as the following: Given an interaction (Hamiltonian), what is the 
most efficient way of entangling particles? Can we make the process more 
efficient by supplementing the action of the Hamiltonian with some local 
unitary operations? Can we increase the entanglement more efficiently by 
using some ancillas? So far, much of the theoretical effort in quantum infor- 
mation theory has been devoted to the characterization and quantification 
of the entanglement of a given state. Very recently, it has been realized 
that there is a parallel notion of the entanglement in the dynamics of a 
system jB]. In [H], the authors consider the situation that one has a given 
unitary transformation and ask, for example, how much state entangle- 
ment is needed to produce it. Here we focus on a different issue: Given an 
interaction (i.e., a Hamiltonian) how can we make the most effective use of 
it [S]? What we propose here is to define and determine the entanglement 
capabilities of physical processes, in particular, of unitary evolutions [ j54] . 
This is a very relevant problem not only from the theoretical point of view 
but also from the experimental one. Of course, this problem is even more 
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difficult tfian tfie one of quantifying the entanglement of states. In any 
case, then we give the first steps in this direction by considering the case 
in which the physical process is acting on two qubits. From our results it 
turns out that (i) it is more efficient to produce entanglement if initially 
one already has some, (ii) The best initial entanglement is universal, i.e., 
independent of the physical process, (iii) One can improve the performance 
of a physical process by complementing it with fast local operations, (iv) 
One can also improve it (in certain cases) by using auxiliary systems, (v) 
All entangling Hamiltonians can simulate each other and are thus qualita- 
tively equivalent; we also provide an upper bound on the time required for 
one Hamiltonian to simulate another. We consider two qubits interacting 
via a nonlocal Hamiltonian H. We want to determine the most efficient 
way in which we can use such an interaction to produce entanglement. We 
will characterize the entanglement of a state of the qubits at a given time 
t, \ip{t)), by some entanglement measure E. In order to quantify the en- 
tanglement production, we define the entanglement rate F at a particular 
time t of the interaction as follows: 



Tit) = '-^. (1.3) 



This quantity depends on \ip) not only through its entanglement E. 
The goal is then to find the conditions which must be satisfied in order to 
obtain a maximal entanglement rate. In particular, we will be interested in 
determining the following: (i) For any initial entanglement E of the two- 
qubit system, what is the state \ip{t)), say l'?/'^), for which the interaction 
produces the maximal rate F^'? (ii) The maximal achievable entanglement 
rate Tmax ( note that the definitions of F and Tmax are not restricted to 
qubits but are also valid for d- level systems, c/ > 0), Tmax = maxETE and 
the state |V^max) for which F = F^^^;. 

In what follows we provide some essential background to appreciate 
the work presented in chapter 2. We review some definitions and known 
results. Let \il)) be a state of the systems A and B. This state can always 
be written using the Schmidt decomposition f5]. 
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= J2^/X^\(|)^)A^h)B (1-4) 

i 

where {|0i)} and {|?7i)} are orthonormal sets of states, and with A^ > 
^- Xi = 1. The entanglement between A and B is defined as 

Eim■.= -Y,>^^^Og,X, (1.5) 

i 

Reference [42] considers maximizing the rate of increase of entanglement 
when a pure state is acted on by e^"^^^, the evolution according to a time- 
independent Hamiltonian H (we set h = 1 ). We refer to this maximal rate 
as the single-shot entanglement capacity. When no ancillas are used, this 
is given by 

r(i*) r E {e-^^^ Ij;)) - E im 

I rr .= max lim — ^^ (l-o) 

\ij)eHABt^O t 

Here the rate of increasing entanglement is optimized over all possible 
pure initial states of Hab without ancillary systems. In fact, the single- 
shot capacity may be higher if ancillary systems A and B, not acted on by 
H, are used. For this reason, we may consider the alternative single-shot 
entanglement capacity 

r)j' := sup lim — ^ (1.7) 

Men - . ^^0 t 



"■AABB 



Note that in Eqs. ( (11.31) , (11.41 )). the limit is the same from both sides 

(1*) / r(i) 



even though it might be the case that T \j* ^ ^ h ^^ general (and similarly 



for ry). For any two-qubit Hamiltonian i7, Ref. [12] shows that it is 
locally equivalent to a canonical f orra 

^ ^iGi (g) (Ji, llx> lly> \llz\ (1-8) 

i=x,y,z 

In terms of this canonical form, the optimal single-shot entanglement 
capacity of any two-qubit interaction without ancillas is given by 
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.(1*) . 



ri^*^= a (^, + ^,) . (1.9) 



a ■-2 max y^x ( 1 - x) log2 « 1.9123. (1.10) 

^ y J- X / 

where the maximum is obtained at xq ~ 0.9168. In addition, F]^ may 



,(1*) 



be strictly larger than T\j when l/i^l > [42]. Reference [43j considers 



the asymptotic entanglement capacity T jj for an arbitrary Hamiltonian 
H . Tjj is defined as the maximum average rate at which entanglement 
can be produced by using many interacting pairs of systems, in parallel or 
sequentially. These systems may be acted on by arbitrary collective local 
operations (attaching or discarding ancillary systems, unitary transforma- 
tions, and measurements). Furthermore, classical communication between 
A and B and possibly mixed initial states are allowed. Reference [13] 



proves that the asymptotic entanglement capacity in this general setting 
turns out to be just the single-shot capacity in Ref. [42], Th = F^ , for all 
iif , so 

E {e-^^' \^)) - E m) 
Th'-= sup lim — ^ (l-ll) 

Note that the definition of the capacity involves a supremum over both 
all possible states and all possible interaction times, but in fact it can be 
expressed as a supremum over states and a limit as t — )• 0, with the limit 
and the supremum taken in either order. Let \il)) be the optimal input 
in Eq. ( (11.41 ) or (11.81 )). When |?/^) is finite dimensional, the entanglement 
capacity can be achieved [121 33] by first inefficiently generating some EPR 
pairs, and repeating the following three steps: (i) transform nE {\'ip)) EPR 
pairs into \ip)®'^ [521 153], (ii) evolve each \ip) according to H for a short 
time 5t^ and (iii) concentrate the entanglement into n [E {\il))) + 5tT h) 
EPR pairs [52]. It shows that F]^ := F]^ for any two-qubit Hamiltonian 
with canonical form 

K := /i^cr^c 0-2; + ii:,CFy ® ay, idx> l~^y>0 (1-12) 



1.5 Entanglement in indistinguishable particle systems 13 

SO that all three entanglement capacities are equal: 

rK = r^K^ = rP (1.13) 

The optimal input is therefore a two-qubit state, and the optimal proto- 
col applies. In particular, for these Hamiltonians, which include the Ising 
interaction cr^ dz and the anisotropic Heisenberg interaction /i^Ca: 'S>(Jx + 
Hxdy <S> (Ty, entanglement can be optimally generated from a two-qubit ini- 
tial state \iIj) without ancillary systems AB. As mentioned above, this 
result is not generic, since ancillas increase the amount of entanglement 
generated by some two-qubit interactions, such as the isotropic Heisenberg 
interaction fi^dx <S> a^ + fJ'xf^y ^ CTy + fi^^^z ^ ^z [42]. 

1.5 Entanglement in indistinguishable particle sys- 
tems 

Understanding and using entangled states of identical and indistinguish- 
able particles [581 ISSl EDI EH ESI ES, EH E5] generates many questions of 
fundamental nature. 

The first of these problems is that of locality and nature of local oper- 
ations [i58j. for the case of distinguishable particles, by locality, we mean 
Einstein locality |ES], which is identically realized by different parts of a 
quantum system being space-like separated. However, the quantum parti- 
cles making up a system are essentially indistinguishable as long as their 
wave functions overlap, that is, they are at short distances from one an- 
other [59] . Such a situation can arise, for example, in a quantum device 
based on quantum dot technology [EZl EHl E2]- Here qubits are realized 
by the spins of the electrons in a system of quantum dots. The overlap 
between the electron wave functions in different dots can be varied by 
controlling parameters like gate voltages or magnetic fields, which change 
the tunneling amplitudes of the electrons from one dot to the other. For 
non-negligible overlaps, the entanglement between the qubits is then inti- 
mately connected to the electron entanglement, which is essentially that of 
indistinguishable fermions. 
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Second problem is to define separable (and hence entangled) states. For 
N-qubits, the state space is the tensor product of state-spaces of individual 
qubits. A N-qubits pure state is separable if it can be expressed as tensor 
product of individual qubit states [5SI ES]. For indistinguishable parti- 
cles all physical states reside in the subspace of the tensor product state 
space, having appropriate symmetry (antisymmetric (symmetric) subspace 
for fermions (bosons)) [321 E2]- This subspace cannot be expressed as the 
tensor product of the state spaces of individual particles. Thus the usual 
definitions of separable states cannot be implemented here, a way out is to 
view N-particle state as separable if its Fock space representation contains 
a single term 

\nun2....) = l[-^(a])'''\0) (1.14) 



[n. 



where ai{a]) annihilates (creates) a particle with individual particle state 
\i/ji). This is often expressed by saying that the state \i/j) has slater number 
1. Note that the separable state of two indistinguishable fermions namely, 
a] |0) is given, in the first quantized version, as 

^ {|10>-|10>} 



V2 

when viewed as a two qubit state, this is a maximally entangle state! 
However, as a two fermion state, this apparent entanglement comes about 
only as a consequence of the anti symmetrization requirement, such an 
apparent "entanglement" (i.e. the corresponding correlations ) cannot be 
used as a resource in a quantum information processing or quantum com- 
munication task. If this was possible, then the results of a local measure- 
ment on a fermion will be affected by the existence of identical fermions in 
the universe, which is not true. 

Additional correlations in many-fermion systems arise if more than one 
terms occur in Eq. (11.141 ) that is, if the Slater number for the state exceeds 
one. i.e., if there is no single-particle basis such that a given state of N indis- 
tinguishable fermions can be represented as an elementary Slater determi- 
nant (i.e., fully antisymmetric combination of N orthogonal single-particle 
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states). These correlations are the analog of quantum entanglement in sep- 
arated systems and are essential for quantum information processing in non 
separated systems. As an example consider a "swap" process exchanging 
the spin states of electrons on coupled quantum dots by gating the tunnel- 
ing amplitude between them [681 ES] . Before the gate is turned on, the two 
electrons in the neighboring quantum dots are in a state represented by 
a simple Slater determinant and can be regarded as distinguishable since 
they are separated by a large energy barrier. When the barrier is lowered, 
more complex correlations between the electrons due to the dynamics arise. 
Interestingly, as shown in Refs. [681 ES]? during such a process the system 
must necessarily enter a highly correlated state that cannot be represented 
by a single Slater determinant. The final state of the gate operation, how- 
ever, is, similarly as the initial one, essentially given by a single Slater 
determinant. Moreover, by adjusting the gating time appropriately one 
can also perform a "square root of a swap" which turns a single Slater de- 
terminant into a "maximally" correlated state in much the same way [73] . 
At the end of such a process the electrons can again be viewed as effectively 
distinguishable, but are in a maximally entangled state in the usual sense 
of distinguishable separated particles. In this sense the highly correlated 
intermediate state can be viewed as a resource for the production of en- 
tangled states. We expect that similar scenarios apply to other schemes 
of quantum information processing that involve cold particles (bosons or 
fermions) interacting at microscopic distances at which the quantum statis- 
tics becomes essential. For instance, it should be of relevance for quantum 
computing models employing ultra cold atoms in optical lattices [71] or 
ultra cold atoms in arrays of optical micro traps [75] . 

A possible way to model and quantify entanglement in the system of 
N identical and indistinguishable fermions is to map the corresponding 
Fock space to an isomorphic A^-qubit space and measure and monitor the 
N-fermion entanglement in terms of that on the mapped A-qubit space. 
This way was suggested and used by Zanardi [58]. We use this approach 
to clarify the above problems and establish a quantitative measure for the 
entanglement in N-fermion systems. We emphasize that our measure can 
deal with multi-partite entanglement in N-fermion systems; an area almost 
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untouched until now. 

Entanglement in fermionic systems has myriads of applications in the 
area of quantum devices and is also expected to play a fundamental role in 
many physical phenomena like quantum phase transitions, quantum Hall 
effect, and so on [701 H]) involving many-body quantum systems. 



1.6 Quantum discord 

Entanglement in a quantum state of a multipartite quantum system is 
a fundamental paradigm to isolate and understand quantum correlation 
imphed by the state. For pure states and a large class of mixed states, en- 
tanglement corresponds to the non-local quantum correlations which break 
Bell inequalities. However, quantum correlations breaking Bell inequalities 
need not account for all quantum correlations in a composite quantum sys- 
tem in a given state. In order to account for the quantum correlation in a 
given state, we must find some means to divide the total correlation into 
a classical part and a purely quantum part. This is particularly important 
for mixed states, since their quantum correlations are many a time hidden 
by their classical correlations (CC). An answer to this requirement is given 
by quantum discord (QD), a measure of the quantumness of correlations 
introduced in Ref. [TB]. Quantum discord is built on the fact that two 
classically equivalent ways of defining the mutual information turn out to 
be inequivalent in the quantum domain. In addition to its conceptual role, 
some recent results [77], suggest that quantum discord and not entangle- 
ment may be responsible for the efficiency of a mixed state based quantum 
computer. 

The idea is to take advantage of the observation that in pursuing quan- 
tum analogs of classical notions, equivalent classical expressions often lead 
to different quantum analogs due to non commutativity of operators which 
represent quantum states and observables, and this difference can be ex- 
ploited to characterize and quantify the "quantumness" of an object. In 
particular, Olliver and Zurek defined quantum discord, as the difference of 
two natural quantum extensions of the classical mutual information, and 
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exhibited its applications in revealing quantum aspect of correlations in 
bipartite states including separable ones. The quantum discord is further 
used by Zurek in analyzing Maxwells demons |78j. A closely related and 
important quantity has also been introduced by Henderson and Vedral 
from a different perspective ||79]. Other similar quantities with the same 
spirit have been extensively studied by Horodecki et al. [80]. 

Correlations between two random variables of classical systems A and B 
are in information theory quantified by the mutual information I{A : B) = 
H{A)+H{B) — H{A, B). If A and B are classical systems, then H{.) stands 
for the Shannon entropy H{p) = —^iPilog2Pi, where p = {pi,P2^---)^ is 
the probability distribution vector, while iif (., .) is the Shannon entropy 
of the joint probability distribution pij. For quantum systems A and S, 
function H{.) denotes the von Neumann entropy S{p) = —Trp log2 p, where 
p is the density matrix. In the classical case, we can use the Bayes rule 
and find an equivalent expression for the mutual information I[A : B) = 
H{A) — H{A\B)^ where H{A\B) is the Shannon entropy of A conditioned 
on the measurement outcome on B. For quantum systems, this quantity 
is different from the first expression for the mutual information and the 
difference defines the quantum discord. Consider a quantum composite 
system defined by the Hilbert space Hab = 'Ha®'Hb- Let dimensions of the 
local Hilbert spaces be dimHA = dA and dirriTiB = dB, while dirnHAB = 
dAB- Given a state p (density matrix) of a composite system, the total 
amount of correlations is quantified by quantum mutual information jHl] : 



I{p) = S{pa) + S{pb)-S{p) (1.15) 

where S {p) is the von Neumann entropy and pa^b = Ttb^a (p) are reduced 
density matrices. A generalization of the classical conditional entropy is 
S (yPB\A)i where Pb\a is the state of B given a measurement on ^.By opti- 
mizing over all possible measurements in A, we define an alternative version 
of the mutual information 



{E,} ^ 
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where pb\a = Tta {Ek Ibp) /Tr {Ek ® Ibp) is the state of B conditioned 
on outcome k in ^4, and Ea represents the set of positive operator valued 
measure elements. The discrepancy between the two measures of informa- 
tion defines the quantum discord [761 [72] : 

Da{p) = I{p)-Q{p) (1.17) 

The discord is always non-negative jTG] and reaches zero for the classi- 
cally correlated states [7^ . Note that discord is not a symmetric quantity 
Da (p) 7^ Db (p) and Da (p) refers to the "left" discord, while Db (p) refers 
to the "righf discord. The state p for which Da (p) = Db (p) = is com- 
pletely classically correlated in the sense of [H21 [Bl] • In this thesis, when we 
refer to the discord we mean the "/e/f discord Da {p)- To give an exam- 
ple of a state with nonvanishing discord, consider the two-qubit separable 
state in which four nonorthogonal states of one qubit are correlated with 
four nonorthogonal states of the second qubit: 

J (|0) (0| ® 1+) (+1 + |1) (1| ® I-) (-1 + 1+) (+1 ® |1) (1| + 1-) (-1 ® |0) (0|) 

(1.18) 
Unlike the state above, one can show that the state {p) is of zero discord 
if and only if there exists a von Neumann measurement 11^ = \ipk) {ipk\ such 
that EM 



E 



/Uk^lB)p{Ilk^lB) (1.19) 

k 

In other words, the zero-discord state is of the form '}2kPk \''Pk) {'^k\^Pk^ 
where \i/jk) is some orthonormal basis set, pk are the quantum states in 5, 
and Pk are nonnegative numbers such that ^^Pk = 1- 

Easily implementable necessary and sufficient condition [83]. Let us 
choose basis sets in local Hilbert-Schmidt spaces of Hermitian operators. 
An and Bn where n = 1, ...,d\ and m = 1, ..., c/^. We decompose the state 
p of the composite system into p = '^^m "^nm^n ® An- The coefficients rnm 
define d\ x d\ real matrix i?, which we call the correlation matrix. We can 
find its singular value decomposition (SVD), URW^ = diag[ci, C2, ...] where 
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U and W are d\ x d\ and (i| x (f^ orthogonal matrices, respectively, while 
diag[ci, C2, ...] is d\xd% diagonal matrix. SVD defines the new basis in local 
spaces Sn = X^^ Unn^n and Fjn = ^rh ^mmBrh- The statc p in the new 
basis is of the form p = ^^^i CnSn ® F^, where L = rankR is the rank of 
correlation matrix R (the number of nonzero eigenvalues c^). The necessary 
and sufficient condition Eq. ( 11.191 ) becomes ^^=1 ^n (Xlfc ^kSrJ^k) ® Fn = 
X^n=i ^nSn ® Fn and it is equivalent to the set of conditions: 

Y,^kSJ^k = Sn^ n = l,...,L, (1.20) 

k 

or equivalently [Sn-, H^] = 0; for all k, n. This means that the set of 
operators Sn has a common eigenbasis defined by the set of projectors 11^. 
Therefore, the set 11^ exists if and only if 

[Sn,S,n] = ^, n, m = l,...,L. (1.21) 

In order to show zero discord we have to check at most L{L — l)/2 
commutators, where L = rankR < min \^d\^ d\^ . Now, recall that the 
state of zero discord is of the form p = ^^LiPk^k ® Pk] therefore, it is 
a sum of at most dA product operators. This bounds the rank of the 
correlation tensor to L < dA- Thus, the rank of the correlation tensor is 
the simple discord witness: If L > dA^ the state has a nonzero discord. 
A correlation matrix can be obtained directly by simple measurements 
usually involved in quantum state tomography. However, the detection of 
nonzero discord does not necessarily require measurement of all (d^c/^)^ 
elements of the correlation matrix (full state tomography). It is sufficient 
that the experimentalist measures that many elements of the correlation 
matrix until he finds d^i + 1 linearly independent rows (or columns) of the 
correlation matrix. 

Geometric measure of discord jSB] . Evaluation of quantum discord given 
by Eq. (11.171 ) in general requires considerable numerical minimization. 
Different measures of quantum discord [87] and their extensions to multi- 
partite systems j84] have been proposed. However, analytical expression 



are known only for certain classes of states [86] . Here we use the following 
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geometric measure, proposed in Ref. [83]. 

Df {p) = mm\\p - x\\\ (1.22) 

where Hq denotes the set of zero-discord states and is the square norm in 
the ||X — y|| = Tr {X — Y) Hilbert-Schmidt space and x denotes the 
classical state. We will show how to evaluate this quantity for an arbitrary 
two-qubit state. 

Two-qubit case [32\'- Consider the case Ha = "H^ = C. We write a state 
p in Bloch representation 

^ / 3 3 3 \ 

/9=-/(g)/ + ^ x^a^ (g) / + ^ ?/,/ (g) cTj + ^ T,ja^ aj (1.23) 

\ i—l 1=1 *J=1 / 

where Xi = Trp{ai<S> I), Vi = Trp{I di) are component of the local 
Bloch vectors, Tij = Trp {ai <S> (Jj) are component of the correlation tensor, 
and (Ji, i G [1, 2, 3], are the three Pauli matrices. To each state p we asso- 
ciate the triple {x, y^ T}. Now, we characterize the set ilo- A zero discord 
state is of the form X = Pi iV^i) {'^i\®pi+P2 |V^2) (^2|<S)/?2, where {|?/^i) , |'02)} 
is a single-qubit orthonormal basis, pi^2 are 2x2 density matrices, and pi^2 
are non-negative numbers such that ^i +^2 = 1- 
We defined = Pi— P2 and three vectors 

e={iJi\d\iJi). (1.24) 

s± = Tr {pipi ± P2P2) ^- (1-25) 

It can easily be shown that te and 5+ represent the local Bloch vec- 
tors of the first and second qubit, respectively, while the vector si is di- 
rectly related to the correlation tensor which is of the product form T = 
es^. Therefore, a state of zero discord x has Bloch representation x ~ 
{te, sV, T = es^}, where \\e\\ = 1, ||s±|| < 1, and t G [—1, 1]. The distance 
between state p and x is given by 
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\\p - X\? = \\pf - '^Trpx + 11x11^ (1-26) 

11 I ll~*l|2 I 11— »||2 I II rri 1 1 2 \ 

= ^ (^1 + IfII + \\y\\ + ll^ll J 

-- (1 + txe + ysX + eTs_) 
+- (l + t^+ ||s+||V||s_| 



4 

where ||r|| = TrT^T. First, we optimize the distance over parameters 
t and Spm- It is straightforward to see that its Hessian is a positive and 
nonsingular matrix. Therefore, the function has a unique global minimum. 
The minimum occurs when the derivative is zero: 

II - 11^ 1 

"^ ^" =-{-xe + t) = 0. (1.27) 



dt 2 

'^"^" U-ye + s+)=0. (1.28) 



ds+ 2 

'^"^" = i (-T^e + s_) = 0. (1.29) 



5s_ 2 

which gives the solution t = xe, s_ = y, s+ = T^e. Since the solution 
lies within the range of parameter, \xe\ , ||^|| , ||r-^e|| < 1, it represents the 
global minimum. After substituting the solution we obtain 



|2 1 ii^i|2 



IP -XI, 4 



X 



+ ||r||^ - e{xx^ + TT^) e] , (1.30) 



which attains the minimum when e is an eigenvector of matrix K = xx^ + 
TT^ for the largest eigenvalue. Therefore, we have 

dT ip) = \ {M\' + ll^ll' - kma.) (1.31) 

where kmax is the largest eigenvalue of matrix K = xx^ + TT^. 
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1.7 Bloch Representation 

Throughout this thesis we use the geometric approach to a density matrix 
via its Bloch representation. The determination of a state on the basis of 
the actual measurement (experimental data) is important both for experi- 
mentalists and theoreticians. In classical physics, it is trivial because there 
is a one-to-one correspondence between the state and the actual measure- 
ment. On the other hand, in quantum mechanics, where a density matrix 
is used to describe the state, it is generally nontrivial to connect them 
P, ES, ESI Sni EI]- the Bloch representation of the density matrix can 
be constructed experimently giving the required connection between the 
density matrix and experiments. 

A/^-level quantum states are described by density operators, i.e. unit 
trace Hermitian positive semidefinite linear operators, which act on the 
Hilbert space Ti ~ C^. The Hermitian operators acting on l-L constitute a 
Hilbert space themselves, the so-called Hilbert-Schmidt space denoted by 
HS{1-L), with inner product {p,(j)hs = Tr{p^a). Accordingly, the density 
operators can be expanded by any basis of this space. In particular, we can 
choose to expand p in terms of the identity operator In and the traceless 
Hermitian generators of SU{N) A^ (i = 1, 2, • • • , A^^ — 1), 



i , 
P 



^(/^+^nA,). (1.32) 

The generators of SU{N) satisfy the orthogonality relation 



{K,Xj)HS = Tr{\Xj) = 26,j, (1.33) 

and they are characterized by the structure constants of the corresponding 
Lie algebra, fijk and gijk, which are, respectively, completely antisymmetric 
and completely symmetric, 

2 

XiXj = —5ijlN + ifijkXk + QijkXk- (1-34) 

The generators can be easily constructed from any orthonormal basis 
{\3)}^=Q in H [i92j. The (orthogonal) generators are given by 



1.7 Bloch Representation 23 

{^*}£r^ = {ujk, Vjk, wi), (1.35) 

when I = 1,- ■ ■ ,N — 1 



T~ ' 



K = wi = Jj^j^Y.^\j){j\ - l\l + l){l + 1|), 1 < / < iV - 1, (1.36) 

while for z = iV, • • • , (iV + 2)(iV - l)/2 

K = Ujk = \j){k\ + |/c)(j|, 
and for i = N{N + l)/2, ■■■ ,N^-1 

K = VJk = -^{\J){k\-\k){J\), 

l<j<k<N. 

The orthogonahty relation Eq. (11.331) implies that the coefficients in Eq. 

(11.321) are given by 

U = —Tr{p\i). 

Notice that the coefficient of In is fixed due to the unit trace condition. 
The vector r = (rir2 • • • Tj^[2_iy G R^ ^^, which completely characterizes 
the density operator, is called Bloch vector or coherence vector. The rep- 
resentation Eq. (11.321 ) was introduced by Bloch [l92j in the N = 2 case and 
generalized to arbitrary dimensions in [22] . Any density matrix in two- level 
systems turns out to be characterized uniquely by a three-dimensional real 
vector where the length satisfies 



\X\ = VXiXi<l. (1.37) 

Therefore, if we define the Bloch- vector space i?(R^) as a ball with 
radius 1: 

5(M3) = {a=(Ai,A2,A3)gM':|A|<1}, 

its element gives an equivalent description of the density matrix with the 
following bijection (one-to-one and onto) map from 5(]R^) to the set of 
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density matrices. 

A — > p = -h + -Xtd, 

B(R^) is called the Bloch ball, its surface the Bloch sphere and its element 
the Bloch vector. The equality in Eq. (11.371 ) (i.e., |A| = 1), the surface of 
the ball (the Bloch sphere) which constitutes the set of extreme points of 
Bloch ball, corresponds to the set of pure states, the points interior to the 
Bloch ball correspond to mixed states. It has an interesting appeal from 
the experimentalist point of view, since in this way it becomes clear how 
the density operator can be constructed from the expectation values of the 
operators Aj, 

{X,)=Tr{pX,) = ^U. (1.38) 

As we have seen, every density operator admits a representation as in 
Eq. (11.321 ); however, the converse is not true. A matrix of the form Eq. 
(11.321 ) is of unit trace and Hermitian, but it might not be positive semidef- 
inite, so to guarantee this property further restrictions must be added to 
the coherence vector. The set of all the Bloch vectors that constitute a 
density operator is known as the Bloch- vector space 5(IR^ ^^). from above 
discussion it is known that in the case N = 2 this space equals the unit ball 
in M.^ and pure states are represented by vectors on the unit sphere. The 
problem of determining i?(IR^ ~^) when A^ > 3 is still open and a subject 
of current research jM]- However, many of its properties are known. For 
instance, for pure states {p^ = p) it must hold 



||r||2 = y ^^^ ^\ urjg,jk = {N-2)rk, (1.39) 

where ||.||2 is the Euclidean norm on M^ ^^. In the case of mixed states, 
the conditions that the coherence vector must satisfy in order to represent a 
density operator have been recently provided in [251 |2S] • Regretfully, their 
mathematical expression is rather cumbersome. It is also known jWl EB] 

that B{R^'-^) is a subset of the baU Dr{R^'-^) of radius R = ^M^^ 
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which is the minimum baU containing it, and that the ball Dr(R ) of 
radius r = \ 2(n~i) ^^ included in BiM.^ ^^). that is, 

Dr{n^"-^) c 5(M^'-i) c Dr{W^"-^). (1.40) 

In the case of bipartite quantum systems of dimensions M x N {% c:=l 
C"^® C^) composed of subsystems A and 5, we can analogously represent 
the density operators as 



tj 



where A^ (A-,) are the generators of SU{M) {SU{N)). Notice that r G 
M^ ^^ and s G M^ ^^ are the coherence vectors of the subsystems, so that 
they can be determined locally. 



PA = TtbP = -j-f{hi + ^ ^A,), pb = TvAp = T^l^iv + ^ s^\). (1.42) 

i i 

The coefficients tij, responsible for the possible correlations, form the 
real matrix T G M(M^ ~ 1) x (^^ ~ 1)? ^^^^ ^is before, they can be easily 
obtained by Uj = ^Tr{p\ ® X^) = ^(A, ® Xj). 

1.8 Chapterwise Summary 

This thesis is concerned with: (a) Production of genuine multipartite entan- 
glement quantum systems (> 2 qubits) as well as bipartite entangled quan- 
tum systems (2 qubits and qutrits). (b) Quantification of both bipartite 
entanglement and genuine multipartite entanglement for indistinguishable 
spin fermions. (c) Thermal quantum discord and classical correlations in 
a two qubit XX model (Heisenberg chain) in a non-uniform external mag- 
netic field, (d) Finding a geometric measure of quantum discord for an 
arbitrary state of a bipartite quantum system. The chapters are arranged 
as follows : 
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Chapter 2 : In this chapter, we develop a geometric approach to quan- 
tify the capability of creating entanglement for a general physical inter- 
action acting on two qubits, two qutrits and three qubits. We use the 
entanglement measure proposed by A. S. Hassan and P. S. Joag for A^- 
qubit pure states (Phys. Rev. A 77, 062334 (2008)). This geometric 
method has the distinct advantage that it gives the experimentally im- 
plementable criteria to ensure the optimal entanglement production rate 
without requiring a detailed knowledge of the state of the two qubit sys- 
tem. For the production of entanglement in practice, we need criteria for 
optimal entanglement production which can be checked in situ without any 
need to know the state, as experimentally finding out the state of a quan- 
tum system is generally a formidable task. Further, we use our method 
to quantify the entanglement capacity in higher level and multipartite sys- 
tems. We quantify the entanglement capacity for two qutrits and find 
the maximal entanglement generation rate and the corresponding state for 
the general isotropic interaction between qutrits, using the entanglement 
measure of A-qudit pure states proposed by A. S. Hassan and P. S. Joag 
(Phys. Rev. A 80, 042302 (2009)). Next we quantify the genuine three 
qubit entanglement capacity for a general interaction between qubits. We 
obtain the maximum entanglement generation rate and the corresponding 
three qubit state for a general isotropic interaction between qubits. The 
state maximizing the entanglement generation rate is of the GHZ class. To 
the best of our knowledge, the entanglement capacities for two qutrit and 
three qubit systems have not been reported earlier. 

Chapter 3 : This chapter reports our work on multipartite entangle- 
ment in a system consisting of indistinguishable fermions. Specifically, we 
have proposed a geometric entanglement measure for N spin- 2 fermions 
distributed over 2L modes (single particle states). The measure is defined 
on the 2L qubit space isomorphic to the Fock space for 2L single particle 
states. This entanglement measure is defined for a given partition of 2L 
modes containing m > 2 subsets. Thus this measure applies to m < 2L 
partite fermionic system where L is any finite number, giving the number of 
sites. The Hilbert spaces associated with these subsets may have different 
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dimensions. Further, we have defined the local quantum operations with 
respect to a given partition of modes. This definition is generic and unifies 
different ways of dividing a fermionic system into subsystems. We have 
shown, using a representative case, that the geometric measure is invari- 
ant under local unitaries corresponding to a given partition. We explicitly 
demonstrate the use of the measure to calculate multipartite entanglement 
in some correlated electron systems. To the best of our knowledge, there is 
no usable entanglement measure of m > 3 partite fermionic systems in the 
literature, so that this is the first measure of multipartite entanglement for 
fermionic systems going beyond the bipartite and tripartite cases. 

Chapter 4 : In this chapter, we investigate how thermal quantum dis- 
cord {QD) and classical correlations {CC) of a two qubit one-dimensional 
XX Heisenberg chain in thermal equilibrium depend on temperature of 
the bath as well as on nonuniform external magnetic fields applied to two 
qubits and varied separately. We show that the behavior of QD differs 
in many unexpected ways from thermal entanglement [EOF). For the 
nonuniform case, {Bi = —B2) we find that QD and CC are equal for all 
values of {Bi = —B2) and for different temperatures. We show that, in 
this case, the thermal states of the system belong to a class of mixed states 
and satisfy certain conditions under which QD and CC are equal. The 
specification of this class and the corresponding conditions are completely 
general and apply to any quantum system in a state in this class and sat- 
isfying these conditions. We further find that the relative contributions of 
QD and CC can be controlled easily by changing the relative magnitudes 
of Bi and B2. Finally, we connect our results with the monogamy relations 
between the EOF, classical correlations and the quantum discord of two 
qubits and the environment. 

Chapter 5 : Quantum discord, as introduced by Olliver and Zurek 
[Phys. Rev. Lett. 88, 017901 (2001)], is a measure of the discrepancy be- 
tween quantum versions of two classically equivalent expressions for mutual 
information. Dakic, Vedral, and Brukner [Phys. Rev. Lett. 105, 190502 
(2010)] introduced a geometric measure of quantum discord and derived 
an explicit formula for any two-qubit state. Luo and Fu [Phys. Rev. A 82, 
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034302 (2010)] introduced another (equivalent) form for geometric measure 
of quantum discord. We find an exact formula for the geometric measure 
of quantum discord for an arbitrary state of a m x n bipartite quantum 
system, using the form for geometric measure of quantum discord given by 
Luo and Wu. 

Chapter 6 : In this chapter, we summarize the work presented in this 
thesis and give the possible ways in which this work may be developed 
further. 



Chapter 2 

Entanglement capacity of nonlocal 
Hamiltonians: A geometric approach 



=4tc1 3T^ cT3r U >f1 I I 

Yoga Sutras of Patanjali 

Creativity of such a person gets one with the universal creation, transcending all limitations of an 

individual. 

In this chapter, We develop a geometric approach to quantify the capa- 
bihty of creating entanglement for a general physical interaction acting on 
two qubits. We use the geometric entanglement measure for N-qubit pure 
states PSj- This geometric method has the distinct advantage that it gives 
the experimentally implementable criteria to ensure the optimal entangle- 
ment production rate without requiring a detailed knowledge of the state of 
the two qubit system. For the production of entanglement in practice, we 
need criteria for optimal entanglement production, which can be checked 
in situ without any need to know the state, as experimentally finding out 
the state of a quantum system is generally a formidable task. Further, 
we use our method to quantify the entanglement capacity in higher level 
and multipartite systems. We quantify the entanglement capacity for two 
qutrits and find the maximal entanglement generation rate and the corre- 
sponding state for the general isotropic interaction between qutrits, using 



the entanglement measure of N-qudit pure states flOOj . Next we quantify 
the genuine three qubit entanglement capacity for a general interaction 
between qubits. We obtain the maximum entanglement generation rate 
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and the corresponding three qubit state for a general isotropic interaction 
between qubits. The state maximizing the entanglement generation rate 
is of the Greenberger-Horne-Zeilinger class. To the best of our knowledge, 
the entanglement capacities for two qutrit and three qubit systems have 
not been reported earlier. 

We have already seen in chap.l that a quantum system evolves to gen- 
erate entanglement provided its parts interact. For such an interaction, 
the Hamiltonian of the total system is not just a sum of the Hamiltonians 
pertaining to each part (local Hamiltonians). Thus, for a bipartite system 
AB, Hab 7^ Ha + Hb but has a term which couples the two parts A and 
B. Together with local operations, this coupling can be used to generate 
entanglement [121 HOlj Hi] , to transmit classical and quantum information 
llOll 1451 11021 1103] and more generally, to simulate the dynamics of some 



other Hamiltonian (say H'ab) and thus to perform arbitrary unitary gates 
on the composite space Hab = 'Ha®'Hb llHl [IMl SS] • 

A lot of experimental work is devoted to creating entangled states of 
quantum systems, including those in quantum optics, nuclear magnetic 
resonance and condensed matter physics [50] . Determining the ability of a 
system to create entangled states provides a benchmark of the "quantum- 
ness" of the system. Furthermore, such states can ultimately be put to 
some information processing task like superdense coding [13] , or quantum 
teleportation [105] . 

The theory of optimal entanglement generation can be approached in 
different ways. Ref. [42j considers single shot capacities. For two qubit 
interaction, without any ancilla qubits, Ref. [42] presents a closed form ex- 
pression for the entangling capability and optimal protocols by which it 
can be achieved. In contrast, Ref. [101] considers the asymptotic entangle- 



ment capacity, allowing the use of ancillary systems and shows that when 
ancillas are allowed, the single shot and asymptotic capacities are in fact 
the same. However, such capacities are difficult to calculate because the 
ancillary systems may be arbitrarily large. In this chapter we exclusively 
deal with the single shot entanglement capacity. Throughout this chapter, 
we take ^ = 1. 

The chapter is organized as follows. In section 2.1 we deal with entan- 
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glement capacity for two qubit states while in section 2.2 we deal with this 
problem involving two qutrits. In section 2.3 we address the problem of 
the entanglement capacity involving the genuine tripartite entanglement 
for three qubits. The discussion of the results for two qubit, two qutrit 
and three qubit cases is included separately in sections 2.1, 2.2 and 2.3 
respectively. 

2.1 The two qubit case 

We develop a geometric approach to calculate the entanglement capacity 
of any two qubit system ( for the case of pure state) interacting via a 
Hamiltonian which is locally equivalent to 

Hi = iiia^ ®(jf + fi2(J2 ^(^2 + l^s(^i '^ ^f • (2-1) 

Here /ii > /i2 > A/-3. In this section we omit qubit and system identifiers 
A,B and ^5. 

We define the single shot entanglement capacity by 

Efe-'^^ip)) -E{\ij)) 

^max ^ ^^^ lim— ^ ^- . (2.2) 

\iP)enit^o t 



The Hamiltonian Hj in Eq. (122]) is given by Eq. (O). ^(|^)) in Eq 



(12.21) stands for the two qubit pure state entanglement measure given by 
us and is shown to have all the essential (as well as many desirable, e.g., 
superadditivity and continuity) properties expected of a good entanglement 
measure [99j. For a A^-qubit pure state 



where HT^^-*!! is the Hilbert-Schmidt (Euclidean) norm of the A^ way array 
T^^-* occurring in the Bloch representation of |^)('?/^| [991 1106] . 

The scenario we address, is as follows [12]. The idea is to supplement 
the interaction Hamiltonian Hj with appropriate local unitary operations 
in such a way that the state of the qubits at any time t is precisely \t/jE{t))^ 
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for which the increase of entanglement is optimal. In order to construct 
such a procedure, we consider the evolution given by Hj to proceed in 
very small time steps 6t. Let us also assume that the qubits are initially 
disentangled. Using local operations, we can always prepare the state 
\i/jo) that is, the product state which most efficiently becomes entangled 
under the action of Hj. After a time step 6t, the state will change and 
its entanglement will increase to 6E. Then, we use (fast) local unitary 
operations to transform the new state of the qubits into the state \i/jse) for 
which r is optimal. Note that this is always possible, since for qubits all 
states with the same value of E", say 6E, are connected by local unitary 
transformations. By proceeding in the same way after every time step, 
and taking the continuous time limit 6t — > 0, we obtain that the state 
of the qubits at time t is always the optimal one, \i/jE{t))- Obviously, in 
an experimental realization, this procedure requires that we can apply the 
appropriate local transformations in times which are short compared to 
the typical time scale thj associated with i7/, th = {emax — e-min)^^^ where 
Cmax and Cjnin are the maximum and minimum eigenvalues of Hj. Note that 
Eq. (12.21 ) defines the entanglement capacity as the maximum achievable 
entanglement rate for a given system with given interactions. We are also 
interested in finding the state \i/Jmax) for which the entanglement rate is 
maximum, (denoted by F"^"^ in Eq. (12.21 )). 

We consider two qubits interacting via the Hamiltonian Hj in Eq. (12.11 ), 
which represents general interaction between two qubits p2] . First we find 
the entanglement rate F given by 



F = lim 



'E(e-'^^\ij)] -E(|^))' 



f- 



Here lip) is given by a general two qubit state in the Bloch representation 
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p= |'0)(^| = - / / + ^ rfcCTfc / + ^ S^/ (g) CT/ + ^ Tfc/CTfc (g) cr^ 
\ k I k,l 

(2.4) 
where ak,i, k,l = 1, 2, 3 are the Pauh operators. We denote by T = [r^j] the 
correlation matrix occurring in the last term of Eq. ( 12.41 ). Tij are defined 
by 

nj = Tr{a, (g) ajp) = {ip\a, (g) (Jj\ip). (2.5) 

rt and si are the components of the Bloch vectors jl06j of the reduced 
density operators pA and pb respectively, given by 



Vk = Tr{akPA) = ii^Wk ® /|^), (a) (2.6) 

si = Tr{aipB) = {ij\l0ai\tjj). (6) 
We define the entanglement of the state \i/j) as 



£;(l^)) = l|r||-i, (2.7) 

where ||T|| = \l}2,ij^i^ij '^^ ^^^ Euclidean norm of T. For two qubits, 
this measure is related to concurrence ||99] and hence to the Von Neumann 
entropy of the reduced density matrix. 

After finding F we maximize it, using a simple geometric argument. It 
is heartening to see that the scenario described above emerges naturally 
out of this geometric method. 

The entanglement rate F is given by (See Eq. (12.31 )) 

dE _d\\r\\ _ 1 

'J 

with Tij given by Eq. (12.51 ). We evaluate fij as follows 
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We now use the equation of motion , 



■^^ = iHr,pl 



where the Hamiltonian Hj is defined via Eq. (12.11 ), to get fl21] 



dr. 



u 



-iTr{a, aj[Hi, p]) = iTr{Hi[a, aj, p]). 



Substituting p from Eq. (12.41) and using the commutation relations fl21] 



we get, using \oi^ o^ = 2i£ijk(7k, {cr^, aj} = 26ij and the expression of Hj in 
Eq. (12.11 ). after a bit of algebra. 



dr, 



u 



dt 



= -2 



k,n l,n 



This gives 



E 



^3 



^ij^ij 



i,k,n j,k,n 



Thus we get, for the entanglement rate F, 



r = 



I mi 



y^ [{r X T,n)n + (S X Tn:)n] ^r 



Here T^n and Tn-. are, respectively, the n th column and row vectors of 
the correlation matrix T = [t^j] . 

The entanglement generation rate F expressed in Eq. (12.81 ) is obtained 
via the temporal evolution of the initial state by the interaction Hamil- 
tonian Hi. This expression for F does not depend on any local unitary 
transformation applied to a qubit. Following the general scenario described 
above, (see the third paragraph of this section), we now lock on to an in- 
stant of time and apply the local unitary transformations to qubits, in 
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order to find the conditions for optimal F and the corresponding two qubit 
state \iIje)- The experimental meaning of this sentence is described as a 
part of the scenario above. In the geometrical approach we have adopted, 
local unitary transformations amount to rotations of vectors in Eq. ( 12.81 ). 
which are the vectors in the Bloch space of individual qubits. We expect 
the entanglement to remain unultered by the local unitaries, which turns 
out to be the case. The entanglement measure in Eq. ( 12.71 ) is not affected 
by local unitaries, as proved in [W]. 



Obviously, F will be maximum if the components of the vector prod- 
ucts occurring in Eq. ( 12.81 ) are replaced by the magnitudes of these vector 
products and the factors in these products are mutually perpendicular. 
Geometrically, this means that the vector products themselves are in the 
directions of the components occurring in Eq. ( 12.81 ) with the other two 
orthogonal components zero. For example, (r x t-i) is along its first com- 
ponent, i.e. along x axis, with its y and z components zero. Thus, in order 
to maximize the first term in Eq. ( 12.81 ). namely. 



E' 



^X T;n)niJ^n = (r X 7^l)l/il + (f X 7^2)2/i2 + (r X 7^3)3/i3, 



we must have vectors (r x 7^1), (r x 7^2) and (r x 7^3) along x^y^z axes 
respectively. This can be done only when one of the vector products is 
zero. Since fii > [12 ^ A^3, we choose (f x 7^3) = 0. Given the vector (r x f^i) 
along the x axis and the vector (r x 7^2) along the y axis, we can choose 
r to be along the z axis and vectors ti and 7^2 along the y and x axes 
respectively. In exactly the same way, maximization of the second term 
in Eq. ( 12.81 ). X]„(s x r„:)„;U-„, makes the vector s along the z axis and 
vectors fi: and r2: along the y and x axes respectively. Writing explicitly 
the components of the vector products in the expression for F (Eq. ( 12.81 )) 
and putting ri^2 = = si^2 we get, 

2 

F = pz^((-r3r2i - S3ri2)/ii + (r3ri2 + S37"2i)At2)- 

Since we are dealing with the two qubit pure states we have ||r|| = ||s|| 
I112J . so that rs = ±53. Choosing 



36 2. Entanglement capacity of nonlocal Hainiltonians: A geometric approach 

ra = -S3 (2.9) 

we get, 

2 

r = -r-^rsiru - r2i)(^i + ^2)- 

The expression (ri2 — T21) becomes maximum when 

ri2 = -rsi. (2.10) 

Finally, we note that this maximization procedure does not change 
ll'^ll ~ V^n \\^-n\\^ ^^^ hence the entanglement value given by Eq. (12.71) . 
Further, choosing the cross products along their components appearing in 
Eq. (12.81) corresponds to the rotations in Bloch space, generating local 
unitaries on the system. Therefore, the maximum of F over the states with 
same entanglement, that is, F^;, is given by 

4 

Fi^ = — — r3ri2(Aii + A^2)- (2.11) 

To get the state \iI)e) corresponding to F£;, we seek the state satisfying 
conditions Eq. (12.91) and Eq. (12.101) . We start with the general state 



I?/;) = '^- • Q Cij\ij) and calculate T12 and r2i. In order to satisfy Eq. (12.101) . 
the state |?/^) should be 

\i^E) = |coi||01) + z|cio||10) ; Icoip + \ciof = 1, (2.12; 



which is the same as {iIje) obtained in Ref [12] if we identify |coi| = ^/p. 



Further, we can write F£; (Eq. (12.111) ) as the product of two factors 

r^ = f{p)hmax 

with 
and 

m = ^. (2.13) 



2.1 The two qubit case 37 

To get f{p) as a function of j9, we calculate rs, tu and ||T|| using the state 
\i^E) (Eq. ( 121^ ) so that 

^^^^ lirii ^i + 8pii-p) ■ 

Fig. (2.1a) depicts this f{p) verses p, while Fig. (2.1b) plots the analogous 
f{p) obtained using Von Neumann entropy as the entanglement measure. 
Note that f{p) and hence r£; vanishes for the maximally entangled state 
(j9 = I) about which it is antisymmetric /(| + a:) = —f{^ — x). We see that, 
as p increases from to |, F^ > makes the entanglement increase, until 
is maximal at j9 = |, after which F < 0, making entanglement decrease to 
zero as p approaches 1. 



Thus we see that, F^; is the product of the function which depends 
only on the state, (via p) and the factor hmax which depends only on the 
interaction strengths fii and fi2, that is, on the interaction Hamiltonian. 
Note that hmax is independent of the entanglement measure. The form 
of f(j)) for the entanglement measure in Eq. (12.71 ) and that for the Von 
Neumann entropy, (Fig. (2.1b)) also turns out to be the same. To get F"^"^ 
we have to find j^o at which f{p) is maximum. To do this, we invoke the 
relation between fvN{p), which is the analog of f{p) in Eq. (12.11 ) obtained 
via Von Neumann entropy of the reduced density operator of the two qubit 
pure state [42] and f{p) in Eq. (Elj) obtained via Edi/j)) in Eq. (EJ}). This 



IS 

fM,) - m (^/f ) . (2.14) 

where Evn is the Von Neumann entropy of the reduced density operator 
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1.5 




Fig. (2.1a) 



TO 1 




Fig. (2.1b) 

Fig. (2.1): (a) f{p) vs. p for entanglement measure in Eq. (12. 7p and (b) fvNip) vs. p for 

Von Neumann entropy of the reduced density matrix (see text). 
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and E is given by Eq. ( 12.71 ). Maximizing the RHS of Eq. (12.141 ) we get 
po ^ 0.0832217 and E"'"^ « 1.9123. The state \ijmax) corresponding to 
Ymax -g ^YiQ state \iIje) with p = pq. 



Erom the definitions of Tj,- and rfc(Eq. (12.51 ) and (12.61 )). these quantities 



are averages of the Eauh operators in a state 1?/^), which can be obtained 
using experimentaUy measured values of the corresponding operators on 
two qubits. Therefore, the geometric method presented here has the ad- 
vantage that the conditions for E^, Eqs. ( fi.PI J and jJlAd) . can be tested 
experimentally giving us an experimental way to check out whether the sys- 
tem has reached the state \iIje)- The value of E^ can also be experimentally 
estimated via Eq. (12.111 ). for given fii + jj,2. Eurther, the function f{p) can 
be estimated experimentally via Eq. (12.11 ) as the system evolves, under 
the given Hamiltonian, toward {i/Jmax)^ or under the local unitaries toward 
\i/je)- These experimental estimations can be carried out without a detailed 
a priori knowledge of the quantum state at any time during the evolution of 
the two qubit system. These facts can be of great advantage in a practical 
implementation of any scheme to entangle two qubits interacting via some 
Hamiltonian or quantum gates jl07j . Eor the production of entanglement 



in practice, we need criteria for optimal entanglement production which 
can be checked in situ without any need to know the state, as experimen- 
tally finding out the state of a quantum system is generally a formidable 
task. We note that, in order to achieve such an experimental determination 
of optimal entanglement production rate using the model in ref [42] 



we 

have to experimentally obtain the values of the Schmidt coefficients of the 
evolving two qubit state, which requires the experimental determination 
of the two qubit state itself. This requires more experimental effort and 
resources (<i + 1 different joint measurements, d = dimension of the joint 
Hilbert space fl08j ) as compared to measuring the quantities in Eqs. (11.91 ) 
and (11.101 ). which are simply the average values of the Eauli operators in 
the state. 
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2.2 The two qutrit case 

The entanglement measure in Eq. (12.71 ) can be generalized to the A^ qu- 
dit pure states which satisfies all the essential (and many desirable, e.g., 
superadditivity and continuity) properties expected of a good entangle- 
ment measure jlOO] . Therefore, we can use it to obtain the entanglement 
generation rates for the multipartite d level systems. Here we find the en- 
tanglement generation rate for two qutrits (labeled A and B) interacting 
via the Hamiltonian 

8 

^/ = E/^^^>^?' (2.15) 

p=l 

where fipS are the interaction strengths satisfying fik ^ fJ'i ^ov k < l^kj = 
1, . . . , 8. Here Xp, p = 1, . . . ,S are the generators of the SU{3) group 
satisfying Tr{XiXj) = 26ij and are characterized by the structure constants 
of the corresponding Lie algebra, fijk and gijk, which are, respectively, 
completely antisymmetric and completely symmetric. 

2 

^i^j ~ o'^ii^s + ifijkXk + gijkXk- (2.16) 

Other useful relations are 

^ifjki = Tr{[X,,Xk]Xi) (2.17) 

ig,kp = Tr{{\,Xk}Xp). (2.18) 

We give here the generators of SU{3) in the |1), |2), |3) basis fl21j to 
be used below. 

Ai = |l)(2| + |2)(l| 

A2 = -z(|l)(2|-|2)(l|) 

A3 = |l)(l|-|2)(2| 

A4=|l)(3| + |3)(l| 

A5 = -z(|l)(3|-|3)(l|) 

A6 = |2)(3| + |3)(2| 

A7 = -z(|2)(3|-|3)(2|) 

A8 = ^(|l)(l| + |2)(2|-2|3)(3|). 
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The action of these generators on the basis states {|1), |2), |3)} is given 
by the foUowing. 

Ai|l) = |2), Ai|2) = |l), Ai|3) = 
A2|l)=^|2), A2|2) = -z|l), A2|3)=0 
A3|l) = |l), A3|2) = -|2), A3|3)=0 
A4|l) = |3), A4|2)=0, A4|3) = |l) 
A5|l)=^|3), Ai|2)=0, Ai|3) = -z|l) 
A6|l)=0, A6|2) = |3), A6|3) = |2) 
A7|l)=0, A7|2)=z|3), A7|3) = -z|2) 
Asll) = ^|1), A8|2) = ^|2), A8|3) = -^|3). 

We use these equations to get the vectors A"^ and A'^ in M^ whose com- 
ponents are the averages Kf = (V^|Ai ® /|?/^) = Tr{\i ® Ip), i = 1, ... ,8 
and Af = {il)\I (g) A^l?/;) = Tr{1 Xip), i = 1, . . . ,8 respectively, where I?/'), 
{p = |?/^)(?/^|) is a two qutrit pure state (see Eq. (12.201 ) and the discussion 
foUowing it). 

The pure state entanglement for two qutrits is given by jlOO) . 



B(|,A» = lir||-3 

where ||T|| = yYlij^i^ij ^^ ^^^ Euclidean norm of T. The general two 
qutrit pure state p has the following Bloch representation. 

p = \^Jj){^Jj\ 



9 
9 

4 



+ T^^^/tAf^Af), 



kj 



Here (A"^'^) = Tr{pA,BX^'^) with A"^'-^ and pa,b (the reduced density op- 
erator) apply to the qutrit A and B respectively, while Tki = {9/4)Tr{Xf^ 
Xf p). The definitions of Tki and F are 

'•J 
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fij = iTr{H[\®\,p\), 

where we have used the Heisenberg equation of motion as in the two qubit 
case. Using Eq. (12.161 (12.171 ). (12.181 ) and the elements of the tensors fijk 
and Qijk in |121] . we get, after some algebra, the following expression for F 

r = -3 (-^\ Y. ^'pf^^p (^^^^^ + ^p^^?) ■ (2.19) 



Expanding the sum in Eq. (12.191 ) and rearranging, we get. 



r = -3 {^^^ Y.^{S)Y,^A{r:, X A^), + (f,: X \%l (2.20) 

where S runs over the triplets 

(1,4, 7), (2, 1, 6), (3, 1, 5), (3, 2, 4), (2, 5, 7), (3, 7, 6), (5, 4, 6), (3, 6, 8), (2, 5, 8) 

and a{S) has values 1, 1/2, 1/2, 1/2, 1/2, 1/2, 1/2, ^3/2, ^3/2 respectively 
for these triplets, r^ and r^: are the vectors in R^ with p E S where S is one 



■p "'^^^ 'p- 

:a,b 



of the above triplets and the index : varies over a given S for fixed ]7. A ' are 
vectors in M^ respectively comprising the components of A"^'^ indexed by 
one of the triplets S. There are in all 54 terms in Eq. (12.201 ) . Unfortunately, 
all these terms are coupled and a simple geometrical procedure to maximize 
r, as in the two qubit case, seems very difficult. However, it is straightfor- 
ward to maximize F numerically over the coefhcients q^, i^j = 0, 1, 2, by 
expressing all the terms in the expression for F (Eq. (12.201 )) as averages 



in the general two qutrit state \iIj) = X^jo-Cijlu)? hi = 0,1,2. We can 
carry out the numerical maximization for the general Hamiltonian in Eq. 
(12.151 ). where the strengths of interaction /i^ have different values. In that 
case, F does not have the simple structure analogous to F = f{p)hmax as 
in the two qubit case. Therefore, we assume isotropic interactions so that 
all interaction strengths are equal to a common value fi. In this case, F has 
a simple form 
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Therefore, we maximize F assuming the isotropic interactions. The result 
is 

^max _ 3,90495^ 

and the corresponding (normahzed) state is given by 

Coo = -0.28317 + ^0.148948; cqi = -0.433055 + ^0.382479 

C02 = -0.117778 + zO.274948; ciq = 0.0625717 - i0.144584 

cii = 0.102783 - ^0.0787094; cu = -0.340939 - z0.324717 

C20 = 0.25066 - zO.167261; C21 = 0.0344755 - zO.244282 

C22 = 0.227159 - ^0.088347. 

After converting this state to the Schmidt canonical form we get the state 
giving the maximal entanglement generating rate for two qutrits, under a 
general isotropic interaction, as 

l^max) = 0.884297|00) + 0.448838|11) + 0.128697|22). 

We find that E{\iljmax)) = 0.677882. This shows that, in order to increase 
the entanglement of a two qutrit system in an optimal way, it is better 
to start with an initially entangled state rather than a product state, at 
least when all the interaction strengths in Hj (Eq. ( 12.151 )) are equal. We 
also note that the optimal entanglement E{\ijjrnax)) is independent of i7/, 
provided, again, that all interaction strengths in Hj (Eq. (12.151) ) are equal. 



2.3 The three qubit case 

We now deal with the problem of entanglement generation capacity for 
three qubits. We emphasize that this is the capacity to generate genuine 
three qubit entanglement and not the bipartite entanglement between any 
two parts of the three qubit system. In this case also, the entanglement 
measure given by Eq. (12.71 ) can be used as this entanglement measure ap- 
plies to A^-qubit pure states and has all the essential (and many desirable. 



44 2. Entanglement capacity of nonlocal Hainiltonians: A geometric approach 

eg superadditivity and continuity) properties expected of a good entangle- 
ment measure ||99]. For the three qubit case, r in Eq. (12.71 ) is the three 
qubit correlation tensor appearing in the Bloch representation of the state. 
Here r is a three way array while for two qubits r was a matrix. The Bloch 
representation of a general three qubit pure state is 

P= I'0)('0| = -{l®I®I + ^riai®I®I 
+ ^ s„/ (g) dn / + ^ q^I ® I ®am 

n m 

In Im 

+ Y^ t^nrn^I ® ^n ®(7m + Y ^'"^"^^ ®^n® Cjjj . (2.21) 



h 



Here r = [rijk] is a three way array while t^"^ = [t ■ ■ ] are matrices. The 
definitions of various symbols in p are as follows. 



n = Tr{afpA)=Tr{af(^I(S)Ip) 
Sn = Tr{(j^pB)=Tr{I®(j^®I p) 
qm = Tr{agpc) = Tr{l0l0agp), 



hn - 


= Tr{af^aSpAB) = Tr{af®a^^l''p) 


• AC 


= Tr{af agpAc) = Tr{af I^ agp) 


,BC 


= Tria^0agpBc) = Tr{I^0a^0aZp): 




Timn = Tr{af0a^0aZp). 



where Pa,b,c and pab,ac,bc are the appropriate reduced density operators. 
We consider the general interaction between qubits which can be reduced 
by the singular value decomposition to the Hamiltonian 

Hi = Hab + Hac + Hbc (2.22) 
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where 



s=l 



s=l 
3 

C 

s 

s=l 
3 



Hac = J2^^fcTf(^l''^a'^. (2.23) 

It is helpful to imagine that the three spins are at the vortices of a triangle. 
If they are arranged on a line, we expect on physical grounds that one 
of the terms can be neglected in comparison with the other two, as it 
gives the next nearest neighbor interaction. For all subsystems we have 
)Ui > A^2 > /^3- Using the definition of the entanglement generation rate F 
in Eq. (12.31 ) and the definition of the entanglement measure in Eq. (12.71 ) 
we get, 

1 ^ 
Using the Heisenberg equation of motion. 



where the Hamiltonian Hj is defined via Eq.s (12.221 ). (12.231 ). we get 



f^jk = iTr{Hi[a^ aj <S) at, p]). 
We now use the commutator identity 

[A0B^C,D^E^F] = ^{[A, D]®[B, E]^[C, F] + [A, D]0{B, E}0{C, F} 
+{A, D} (g) [B, E] {C, F} + {A, D} {B, E} [C, F]) 
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and the definitions of p and Hj in Eq.s (12.211 ) and ( 12.221 ) respectively to 
get, 

3 3 



'^ijk — ^ 



j'=l fc'=l 

3 3 

^nBC\^ fAC^ ^ ,,AC\^.BC^ , 



^'=1 
3 



3 



'U J 



(2.24) 



fc'=l j'=l 

where e s are the Levi-Civita symbols. Substitution of Eq. (12.241 ) in the 
expression for V gives, 

3 



r = 



I mi 



Y, [{r:sk X tf), + (f,:fc X t"5^) J pf + 



fc,s=l 
3 



+ E [(^" ^ *?'')» + (^'^ X *?'')»] '^?'' + 



i,s=l 
3 



;2.25) 



where T-_sk = [^iski ^2sk, ^Ssk]'^, for example, is a vector in R^ for fixed s and 
k. Similarly, t'.'j^ and t'j. are the /cth column and the jth row vectors of the 
matrix t^"' . The expression for the entanglement generation rate F for three 
qubits (Eq. (12.251 )) has 54 coupled terms, each term being a component of 
the cross product of two vectors. A geometric argument to maximize F, 
as in the two qubit case, seems to be very difficult. However, it is quite 
straightforward to maximize F numerically, by writing the elements of the 
three way array T and the matrices t^ "-' as the appropriate averages in the 
general three qubit state 



= ^q|z) ; ^|qP = 1 

i=0 i 



;2.26) 
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where i labeling the product basis ket \i) is the binary representation of the 
index i. We can numerically optimize F for the general Hamiltonian in Eq. 
(12.221) . However, for the general case, where the interaction is anisotropic, 
that is, the strengths of interaction /i^"^ have different values, F does not 
have the simple structure F = f{p)hmax as in the two qubit case. Therefore, 
we assume isotropic interactions so that all interaction strengths are equal 
to a common value /i. In this case, after evaluating all the terms in Eq. 
(12.251) in the state li/j) given by Eq. (12.261) . F can be written as 

F = /i(co, . . . , 07)11. 

After the numerical optimization of F as a function of q, z = 0, 1, . . . , 7, 
we get, 

j.max ^ 5 J2523^. 

The (normalized) state corresponding to this F"^"^ is given by 

|^^„^) = (0.033768 - ^0.168758|000)) + (0.574022 - z0.0709471)|001) + 

+(0.0218412 - z0.111565)|010) + (0.672021 - z0.0754116)|011) + 

+(-0.0603488 + ^0.172566)|100) + (-0.0051137-^0.183831)1101) + 

+ (0.0556843 + i0.151888)|110) + (0.0700719 - i0.259423)|lll). 

This state has the following Acin canonical form, expressed by the two fold 
degenerate sets of entanglement parameters |109j . 



|^+) = 0.610291|000) + 0.67402 exp(i2.51395)|100) + 0.394893|101) + 

+0.110357|110) + 0.0715772|111), 
or, 

1^-) = 0.329873|000) + 0.546087 exp(^0.402558)| 100) + 0.730583|101) + 

+0.20417|110) + 0.132424|111). 

We see that the state with the maximal entanglement generation rate F 
belongs to the GHZ class. Further, we find that E{\i/jmax)) = 0.258918. 
This means that, given the isotropic interaction, it is beneficial to start 



max 
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with an entangled three qubit state for optimal entanglement generation. 
Also, we note that the optimal entanglement is independent of Hj^ provided 
the corresponding interaction is isotropic. 

Thus we see that, for three qubits, the geometric method based on the 
entanglement measure given by Eq. (12.71 ) can be numerically implemented 
to get the state with maximal entanglement generation rate. This pro- 
gram can be carried out for the general interaction Hamiltonian in Eq. 
( 12.221 ) although we have restricted to the isotropic interactions. This pro- 
cedure can be suitably carried out in a laboratory using quantum circuits. 
Every quantum circuit acts unitarily on a quantum state and we can al- 
ways find a Hamiltonian corresponding to such a circuit [IHl I109J . On the 



other hand, given a (interaction) Hamiltonian for a three qubit system, we 
may construct a circuit implementing the corresponding evolution using 
universal quantum gates. We note that, for the isotropic interaction, the 
maximal entanglement generation rate F"^"^ is proportional to the inter- 
action strength ^ and the corresponding state is independent of /i, as in 
the two qubit case. When the interactions are anisotropic, the scenario for 
two qubits does not apply to the three qubit case, as F does not factor 
into the product of a state dependent function and an expression involving 
only the interaction strengths. Even when the interactions are isotropic, 
we do not know the explicit form of such a state dependent function. In 
other words, we do not know whether it is possible to separately account 
for the contribution due to the state and that due to the interactions. Thus 
a general procedure for the maximization of the entanglement generation 
rate for the higher dimensional and multipartite systems still seems to be 
an open question. These observations ensue from the fact that the terms 
in the expression for F could not be decoupled. This difficulty seems to 
be generic, as it may be a consequence of the difficulties in the geometric 
interpretation of the Bloch space for the multipartite and higher dimen- 
sional systems [ jllOl llllj . All the remarks in this paragraph apply to the 



two qutrit case as well. 



Chapter 3 

Multipartite entanglement in 
fermionic systems via a geometric 
measure 



Before I came here I was confused about this subject. 
Having listened to your lecture I am still confused. But on a higher level. 

Enrico Fermi 



3.1 Introduction 

In this chapter, We study multipartite entanglement in a system consisting 
of indistinguishable fermions. Specifically, we have proposed a geometric 
entanglement measure for A^ spin-| fermions distributed over 2L modes 
(single particle states). The measure is defined on the 2L qubit space iso- 
morphic to the Fock space for 2L single particle states. This entanglement 
measure is defined for a given partition of 2L modes containing m > 2 sub- 
sets. Thus this measure applies to m < 2L partite fermionic system where 
L is any finite number, giving the number of sites. The Hilbert spaces 
associated with these subsets may have different dimensions. Further, we 
have defined the local quantum operations with respect to a given partition 
of modes. This definition is generic and unifies different ways of dividing a 
fermionic system into subsystems. We have shown, using a representative 
case, that the geometric measure is invariant under local unitaries corre- 
sponding to a given partition. We explicitly demonstrate the use of the 
measure to calculate multipartite entanglement in some correlated electron 



50 3. Multipartite entanglement in fermionic systems via a geometric measure 

systems. 

We study multipartite entanglement in a system consisting of identical 
particles. We use the idea due to Zanardi [3B] whereby the Fock space of a 
system of fermions is mapped to the isomorphic qubit or 'mode' space. We 
then discuss entanglement in this 'mode' space via a geometric measure. 
The idea is to use the Bloch representation of the state of the m-partite 
quantum system [94] . The measure is defined by the Euclidean norm of the 
m-partite correlation tensor in the Bloch representation. This correlation 
tensor contains all information of genuine m-partite entanglement (see sec- 
tion [32]) • Such a measure was proposed earlier by P. S. Joag and A. S. M. 
Hassan for A^ qubit and N qudit pure states [SS, llOOj and shown to satisfy 
most of the properties expected of a good measure. The Bloch represen- 
tation of a quantum state has a natural geometric interpretation, which 
is why we call this measure a geometric measure [El]. Other geometric 
measures are based on the distance of the given state from the set of sepa- 
rable states in the Hilbert space [1141 1115] . An important question in the 



context of quantum entanglement is that of locality. For indistinguishable 
particles distributed over 'modes' (which are taken to be single particle 
states of particles constituting the system), local operations have meaning 
only in the context of partitions over modes. A quantum operation con- 
fined to a single subset in a given partition is then a local operation. We 
therefore define entanglement in such a system with respect to partitions 
and require it to be invariant under local unitaries defined with respect 
to a given partition. We explicitly demonstrate the use of the measure to 
calculate entanglement in some correlated electron systems. 

The outline of the chapter is as follows: In Sec J3.2l we begin by briefly 
reviewing some details about the isomorphism between the Fock space of 
a system of indistinguishable particles and the 'mode' space [58]. This will 
help us in defining various quantities and also set up notation necessary 
for the subsequent analysis. In Sec J3.3l we define and construct the geo- 
metric measure. Various properties of the measure are also discussed with 
reference to a specific example in Sec J3.4[ In Sec J3.5l we use the measure 
to study entanglement in the Hubbard dimer and trimer. 
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3.2 Mapping between Fock space and qudit space 

We deal with A^ spin-| fermions on a L site lattice. The total number 
of available (localized) single particle states are then 2L in number. The 
fermionic Fock space in the occupation number representation has basis 
states of the form \n1n2 . . . n2L) (n^ = 0, 1 ; i = 1, . . . , 2L). We further 
assume that the total number of particles is conserved. This means that 
we only deal with subspaces of the Fock space corresponding to a fixed 
eigenvalue for the total number operator. We shall refer to this number 
super-selection rule as N-SSR. For a A^-fermion system, we call such a 
subspace of the Fock space W-sector' and denote it by F^. The A^-sector 
of a 2L mode system is the subspace F/v of dimension (^j of the Fock 
space with the dimension of the Fock space for 2L single particle states 
(iV = 0, 1,2, ...,2L) being 

E Ci) - 2- (3.1) 

Since a 2L qubit Hilbert space (C^)*^^^ has exactly this dimension, it is pos- 
sible to construct an isomorphism between the Fock space and the 2L qubit 
Hilbert space (C^)®^'^ jll6] . The particular isomorphism we implement is 



\n\n2 . . . n^h) -^ 1^1) ® 1^2) ®- ■■® \n2L) ; nj = 0, 1 ; z = 1, 2, . . . , 2L (3.2) 

where, in qubit space we associate |0) o | t) and |1) ^ | \). Note that 
the Slater rank of the Fock basis states \n\n2 ■ ■ ■ tl^l) is 1 so that these are 
separable states. Thus the above isomorphism maps separable basis states 
in Fock space to the separable basis states in qubit space. Further, the 
subspace structure of the Fock space namely, 

2L 
F2L = 0FAr (3.3) 

is carried over to the qubit space under mapping (Eq. (15.21) ) because each 
subspace of the Fock space with conserved fermion number A^ is mapped 
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to a subspace of the 2L qubit space spanned by the basis vectors with A^ 
ones and 2L — N zeros. We can write 

H2L = {Cy^ = ^ H2l{N) (3.4) 



where H2l{N) is the image of F/v in Eq. (15.31 ) under the map given by Eq. 

(E2D. 

Next crucial step is to transfer the action of the creation and annihilation 
operators on Fock space to the qubit space, under the isomorphism given 
by Eq. (15.21 ) [ IB5] . We need the creation and annihilation operators a and 
o) acting on a single qubit state, 



a|0) =0, a\l) = |0) 
a^O) = |1), a^l) =0 



(3.5) 



such that. 



Gi - 


-> /® /® ••• (g) ^_a^ 0---(g)/ 




ith qubit 


4 - 


-> /(g) • •• (g) oj^ g)-- • g)/ 




ith qubit 



(3.6) 

Here a^ {a]) is the annihilation (creation) operator acting on Fock space 
F2L: annihilating (creating) a fermion in ith mode. / is the identity on 
single qubit space. The tensor product satisfying the correspondence in 
Eq. (15.71 ) must be consistent with the anti-commutation property of the 
Fock space creation and annihilation operators, 

{tti, aj} = 6tj {a^, ttj} = = {aj, aj} (3.7) 

This requirement leads to the following action of the tensor product oper- 
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ators on the 2L qubit states 

(/ / • • • (g) a{a^) • • • (g) /)(|ni) g) • • • In^) • • • \n2L)) 

ith place 

= (-l)^^'-+i''-'"(|ni) (g) • • • (g) a{a^)\n,) g) • • • g) \n2L)) 



ith qubit 



(3.8) 



Here Ui G {0, 1} ; i G {1, 2, ... , 2L} and Ylij^i+i^i i^ evaluated mod 2. 
Using Eq. ( 15.101 ) , it is straightforward to see that 



{I®!®- 


•• g) a(a^) g)- 


••g)/ 


, I®I® 


■ ■ g) a(a^) g) • 

jth place 


••g)i 


( 


ni) g) • • • g) n^) 


g)--- 


g) nj) g) • 


■ <S) n2L)) = 




and 












{/g)/g)-- 


• g) a_^ g)- • • 

iift, place 


g)/, 


I ®I ®-- 


ji/i place 


®/} 


( ni) g) • • 


• g) n^) g) • • • g) 


nj) g) • • • g) 7121 


^)) = (;/g)...g) 


_^% 










2i factors 





(3.9) 



(3.10) 



We note that the phase factors appearing in Eq. (15.101 ) are the consequence 
of the conservation of the parity operator jll7] 



i2L 



P = Ufti{l-2a]ai). (3.11) 

Henceforth, in this chapter, by 'fermions' we mean spin-| fermions. Fur- 
ther, we call a single particle state a mode. Thus two spin-| fermions on 
two sites is a four mode system. In general, A^ spin-| fermions on L sites 
is equivalent to A^ fermions on K = 2L modes. For example, two spin-| 
fermions on a two site lattice. A, B say, generate four single particle states 
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or modes 1^4 t), |^ i), 1^ t)? |^ i)- In this work, we deal with entangle- 
ment between subsets forming a partition of a 2L-mode fermionic system. 
We define the entanglement measure for any such partition of a 2L-mode 
system without any restriction on the number and the size of the subsets 
forming the partition. These subsystems may involve different degrees of 
freedom, for example, we can deal with the entanglement between spins 
and sites or entanglement between two spins on the same site (intrasite en- 
tanglement). Or if each of the subsets partitioning the 2L modes comprises 
modes with common site label we have the entanglement between sites or 
the 'site entanglement'. Thus all physically realizable subsystems of A^ 
fermions over 2L single particle states can be addressed by dividing the 
2L modes into suitable partitions, for example the 'particle entanglement' 
defined in I60l. 



We now define the local and non-local operations on the 2L mode 
fermionic system [581 ESI EO, 11181 I119J . We do this by using the corre- 
sponding operations on the isomorphic qubit space H2l- We note that, due 
to isomorphism between F2L (Eq. ( 15.21 )) and H2L (Eq. ( 15.51 )). partitioning 
2L modes is equivalent to the corresponding partitioning of the 2L qubit 
system into subsystems. Locality is defined with respect to the partition 
of 2L qubits (or, the corresponding partition of 2L modes) between whose 
subsets we are seeking entanglement. The operations on the state space 
of a single subset in a partition of 2L qubits is taken to be local. The en- 
tanglement measure defined with respect to a partition must be invariant 
under a unitary operation which is local with respect to that partition. We 
will illustrate this point later, using the geometric entanglement measure 
defined below. Henceforth 'mode' and 'qubit' are taken to be synonymous 
and we shall use the expression 'modes' instead of 'qubits'. In other words, 
the spaces F2L and H2L are taken to be the same. 
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3.3 Geometric measure for entanglement. 

3.3.1 Definition 

We define a geometric measure for ttie partitions of the 2L mode A^ fermion 
systems in pure states. Although the definition of the measure is quite 
general, the fermion number super-selection rule restricts the pure states 
to the appropriate subspace corresponding to A^ fermions, namely H2l{N) 
(Eq. (ES])). Thus a state |0110) G H^{2) may be partitioned as |01) (g) |10) 
or as |011) ® |0) etc where the definition of the tensor product is consistent 
with Eq. ( 15.101 ). We use the Bloch representation of A^— partite states 
(drawn from H2l{N)) to get this measure [Ml llOOj - 

First we assume that a partition equally divides 2L modes into subsets, 
i.e. all subsets in the partition contain equal number of modes, say n. 
This corresponds to the case of H2L divided into subspaces of dimension 
(i = 2", n being some divisor of the number of modes 2L. To get the 
entanglement measure we expand the state p = |'?/')('?/'| of the system, 
supported in the appropriate H2l{N), in its Bloch representation. 

In order to give the Bloch representation of a density operator acting 
on the Hilbert space C^ C'^ • • • C'^ of an m = 2L/n-qudit quantum 
system, we introduce following notation [lOOj . We use k, ki {i = 1,2, . . .) 



to denote a qudit chosen from m qudits, so that k, ki {i = 1,2, . . .) take 
values in the set A/" = {1,2,..., m}. The variables at or at^ for a given k 
or ki span the set of generators of SU{d) group for the kth or kith qudit, 
namely the set {Ai,A2,--- ,Arf2_i} for the kith qudit. For two qudits ki 
and k2 we define 

A£i^) = (/rf (g) /d (g) • • • Xa,^ /d • • • /d) 

A£f = (/d /d • • • Xa,^ 0ld^---^ld) 

^all^a'l = {Id®Id®---® A,,^ ®Id®---® A«,^ ®Id® Id) 

where Aq,^ and Aq,^. occur at the kith and A;2th places (corresponding to 
kith and k2th qudits respectively) in the tensor product and are the ak^th 
and ak^th generators of SU{d), ak^ ^ = 1,2, . . . ,d? — 1 . Then we can write 



56 3. Multipartite entanglement in fermionic systems via a geometric measure 



p = ^{'T + E E «-^£' + E E i«.»^.^i!^i +■■■+ 

{ki,k2,--- ,kM} afeiafe2'""fcM 

+ / ^ '^aia2---aN'^ai'^a2 ' ' ' "^aN S (o.lzj 

aia2---aAf 

where s'^'^^ is a Bloch vector corresponding to kth qudit, s^^^ = [safc]^^^! 
which is a tensor of order one defined by 



where pk is the reduced density matrix jl20] for the kth qudit. Here 



{ki, /c2, • • • , ku}, 2 < M < m, 

is a subset of J\f and can be chosen in (]^*) ways, contributing (]^) terms 
in the sum X^r^^ ^^ ••• /ca/I ^^ ^^' P-121 ). each containing a tensor of order 
M. The total number of terms in the Bloch representation of p is 2"\ We 
denote the tensors occurring in the sum X^r^^ ^^ ••• fcA/j' (2 < M < m) by 

'y-{fci,fc2,-,fcAf} = [t„,^c.fe2-afcAJ ^^^^^ ^^*^ defined by 

where Pkik2...kM i^ the reduced density matrix for the subsystem {/ci, /c2, . . . 
, /cm}- Each of the (]^) tensors of order M, occurring in the Bloch represen- 
tation of /?, contains all information about entanglement of the correspond- 
ing set of M subsystems. All information on the entanglement contained 
in p is coded in the tensors occurring in the Bloch representation of p. The 
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tensor in last term in Eq. ( 13.121 ). we call it T''^\ contains all the informa- 



tion of genuine m-partite entanglement. This follows from the observation 
that all other terms in the Bloch representation of p (Eq. (13.121) ) corre- 
spond to subsystems comprising M < m qudits and the density operator 
contains all possible information about the state of the system. 

The operators A^^,, a^ = 1,2, . . . ,(f — 1 are given by ^121j 

A = {ilu, Wl3, W23, • • • , '?^12, Vi3, V23, . . . ,Wi,Wi, . . . , Wd-l} (3.13) 



with 



Ujk jk I -'fcj 

Vjk = -i{Pjk- Pkj) 



wi = WyTT^:^ (^11 + --- + Pii- iPi+iHi) 



l<j<k<d; 1<1 <d-l 



(3.14) 



where 



Pki = \k){l\ {k,l = 1,2,..., d). 

Note that each of the generators of the SU{d) group A^ , i = 1,2, . . . ,d'^ — l 
acts on a single qudit space and hence is local 13.21 apart from the phase 
factor contributed by their action, as given by Eq. (15.101) . We assume 
these phase factors to be absorbed in the coefficients in the expansion of 
the density operator p. 

Let a 2L mode A^ fermion system be partitioned by m = 2L/n subsets, 
each containing n modes. Then for this partition, we define the entangle- 
ment measure for a state |?/^) G H2l{N) by flOOl 



^= Ikll - Ikllsep (3.15' 

where 






d2-l 



\ E *i«... (3.16) 
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and ||r||sep is ||t|| for separable (product) m qudit state 

Irllsep= ;; (3.17) 



3.3.2 Entanglement in partitions with unequal subsets 

We can also generalize the definition of the entanglement to the case where 
the corresponding qubit subsystems have unequal dimensions. We discuss 
the simplest case of bi-partite entanglement with partitions having unequal 
dimensions, say di and d2. In this case, the definition of the geometric 
entanglement measure generalizes to 

E=\\r\\-\\T\\,ep (3.18) 

where 






dj-idl-i 



;EE*?,- (3-19) 

\ 1=1 J=l 



where 



*" = (1) (t) <^l^.- ® ^^1^') = (1) (y) ^''- (3-20) 

Here \i {i = 1, . . . ,d\ — 1) and Xj (j = 1, . . . , 0^2 — 1) are the generators of 
SU{di) and SU{d2) respectively. ||r5ep|| is given by 

iir i|2-iis(i)|iiis(2)|i- r^i^^^^V^^^^^^^ n2n 



Here ||s*^^^|| and ||s^^)|| are the norms of the Bloch vectors of the reduced 
density operators for each subsystem. 

It is straightforward to extend these definitions to partitions containing 
more than two subsets. 
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3.4 Entanglement in a four mode system 

With the entanglement measure defined as above, we give an example 
wherein we can compute the entanglement for different partitions and il- 
lustrate our comments on local and non-local operations. We also compute 
the upper bounds on the entanglement. 

3.4.1 Local and Non-local operations 

Consider a four mode system and the normalized state |?/^) G H4(2) defined 
as 

\iIj) = ^{za|1100) + |1001) + |0110) + |0011)+;5|0101) + |1010)} a'^+(3^ = 2 

(3.22) 
where a, (5 are real. Note that l?/^) can be treated as a member of the Fock 
space -^4(2) with the kets appearing in it being its basis states. Consider 
the evolution of the system in state 1?/^) G F/^{2) via the Hamiltonian 

H = /(a|a4 + oJ^ai) + qh\h2 + Th\ + 7713 + rj{a\a2 + a\ai) (3.23) 

acting on F4. Here / term is the interaction between two modes on different 
sites (intersite interaction), 7/ term is the interaction between two modes 
on the same site (intrasite interaction). F and 7 correspond to single 
mode on site A and B respectively, q term involves number operators 
fii = a\ai ; z = 1,2 for first two modes, on A site. We have included all 
the different kinds of typical interactions encountered in condensed matter 
systems, respecting number super-selection rule. After an infinitesimal 
unitary evolution via this Hamiltonian, the state \il)) evolves to 

IV^') = 1^) - ze//|^> (3.24) 

By employing the mapping of annihilation and creation operators in Eq. 
(15.71 ) and Eq. (15.101 ) and that of Fock space basis states in Eq. (15.21 ). we 
get, for \il)') 

\iIj') = ^{{ia + ief + age)|1100) + (1 - iTe - ier]^)\1001) 
V6 
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+(1 - i-fc - ier])\OllO) + (1 - lef - ^e7)|0011) 

+ {13- efP - ier])\0101) + (1 + lef - leT - le-f - i€r])\1010)} (3.25) 

Now we find the entanglement for different partitions of this four mode 
system, using the geometric entanglement measure, Eq. ( 15.171 ). We first 
partition four modes into four subsets, each containing one mode. This 
case gives genuine entanglement between four modes, which is more general 
than only the bipartite entanglement considered in the literature. For this 
case d = 2, so that ||r||sep = 1 and we get, for the genuine four mode 
entanglement, 

E=\\r\\-1 (3.26) 

where 






, E *|h (3-27) 

\ i,j,k,l 



d=l 



with 

tijki = Tr[p a^ ® aj (ifc ai] = {i/j\ai aj <^ak<^ cri\i^)- (3.28) 

where {ai} i = 0,1,2,3 are the generators of the SU{2) group (Pauli 
operators). The resulting entanglement in \^JJ') is 

EM)) = 
1 

6 

4 {4fa - 2qa{l + /3) + fa^{a^ - 13^) + 4ar/(l + f3)) e ^ ^^ ^ 



6 + ^88 + 64a2 + 32^ + 64/^2 + lOa^^^ + ^4 

+ 0[e'] (3.29) 



-6 + v^88 + 64a2 + 32/5 + 54^2 + 10^2^2 + ^4 

where the first term gives the entanglement E{\ip)) for the state ^JJ as 
defined in Eq. (15.201 ). For this partition, the operations on a single mode 
are the only local operations, while all others are non-local. Therefore, the 
terms Thi and 777,3 are the only local interactions. Therefore, we expect 
that the four mode genuine entanglement should not depend on F or 7 to 
the first order in e, which is the case, as seen from Eq. (15.221 ). 



3.4 Entanglement in a four mode system 6i 

Next, we consider the partition consisting of two subsets, each contain- 
ing two modes on each site, {^4 t, ^ i} and {B t, ^ 1} (site partition). 
Thus we have two subsystems with d = 4 corresponding to a SU{4)<S)SU{4:) 
qudit system. Further, n = 2 giving m = (2L/n) = 2 so that the geometric 
entanglement is 

E.s{\i^')) = \\r\\ - Q (3.30) 



where 



Irll =4 



\ 



15 



2 

k 



with 



Kjk = (^|Aj AfclV^) 

where Aj ; j = 1, . . . , 15 are the generators of SU{4:). The entanglement of 
\iIj') in Eq. (15.221 ) is then given by 



EsiW)) = 

I (-18 + v^208 + 136^2 + 9^4 _ 32^ + 104^2 + S4a^p^ + 9^4^ 

16{-fa + fa^{a^-^'-2))e _^ , 

3(^208 + 136^2 + 9^4 _ 32^ + 104^2 + Ma^/S^ + 9/^4) 

According to the 'site partition', in addition to the operations on single 
modes, the operations on the pair of modes having the same site label are 
also local. Therefore, the resulting entanglement cannot change under the 
intrasite operations in the Hamiltonian, namely the q term, the 77 term and 
as before, F and 7 terms. Thus, to the first order in e, the entanglement 
is expected to depend only on the non-local part of the Hamiltonian, that 
is, on the / parameter. From Eq. ( 15.291 ) we see that this is the case. The 
intersite entanglement quantified using Von-Neumann entropy has been 
reported earlier [58]. 

Thus we see that the geometric measure has the capability to quan- 
tify the genuine multi-mode fermionic entanglement as against the mainly 
bipartite entanglement reported in the literature. Also, the geometric en- 
tanglement measure, for the given partition of modes, shows the correct 
behavior under local and non-local unitary operations. 
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3.4.2 Upper bound for the geometric measure in a four mode 
system. 

We discuss here some upper bounds that one can find for the entanglement 
for the four mode system and compare with existing results obtained from 
other measures. We find that the general state which leads to a maximum 
intersite entanglement computed via the geometric measure also leads to 
a maximum for the von-Neumann entropy. 

We treat the entanglement for the given partition to be the function of 
the coefficients of the general state \il)) G H^{2), namely, 

l^)= E ^^1^) (3.32) 

A:=3,5,6,9,10,12 

where each ket is labeled by the (four bit) binary representation of k. We 
then maximize E{\il))) with respect to the coefficients c^. 

For the site partition (see above) the entanglement (Eq. (15.281) and two 
equations following it) is a function of the coefficients Ck in state \iIj) G 
H/^{2) given in Eq. (13.321) . We find that the the maximum value of the 
entanglement is given by 

E^a. = 1.74593 (3.33) 

The Schmidt canonical form (with respect to the H^{2) basis) of \'ilJmax) 
corresponding to this Ej^ax is 

\i^^ax) = \[ Y. 1^)) (3.34) 

\/i:=3,5,10,12 J 

or, 

\fc=3,6,9,12 J 

where each ket is labeled by the four bit binary representation of k. Also, 
both these forms lead to the von-Neumann entropy = 2 which is the max- 
imum possible 1122] • 
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For the four mode entanglement (Eq. (15.231 )) we find that entanglement 
is maximum for the state given by 

C3 = 0.06701 - 0.33751^; C5 = 0.16442 - 0.51281^ 

cq = -0.1006 + 0.2854^; cg = 0.29967 - 0.04206z 

cio = -0.53522 + 0.05955z; cu = -0.34410 - 0.00118z 

(3.36) 

and the maximum value is found to be 

However, we do not have any entanglement measure to compare with the 
geometric measure. We also note that, for four modes, no canonical form 
like the Schmidt or Acin canonical form (for two and three modes respec- 
tively) is available. 

3.5 Applications 

We now consider some correlated lattice models and discuss multi-mode 
entanglement in these models using the geometric measure. 

3.5.1 Hubbard dimer. 

The Hubbard dimer model is a simple model for a number of physical 
systems, including the electrons in a H2 molecule, double quantum dots, 
etc [581 EQl 1123] . The Hamiltonian can be written as 

o-=t,i 3=A,B 

where ^4, B are the site labels and t and 4- are spin labels, t is the hopping 
coefficient measuring hopping between two sites while conserving spin and 
U quantifies Coulomb interaction between fermions on the same site. By 
varying (^) we can vary the relative contributions of hopping and Coulomb 
mechanisms. 
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The ground state of the system at zero temperature can be easily ob- 
tained as [581 ED] 

l^o) = NGQ\vac) (3.39) 

where N = (?/^o|'^o)~^ is the normahzation factor and 



.t .t 



.t J 



'U' 



Go — Cj^^Cj^^^ + c^^c^^ "^ *^ I 4^ ) y^A-f^Bi ^Ai^B-t 



J J 



J J 



(3.40) 



with a{x) = X + vT+x^. By mapping to H2L via Eq. (15.21 ) we get, 

\vac) -^ |0) |0) |0) (g) |0) = |0000) (3.4i; 

while mapping between the operators, (Eq. (15.71 ) and Eq. (15.101 )) gives 



--At 



o) ®h®h®h c 



,t 

Ai 



h^a^ ®l2® h 



l2®h®o)®h c^Bi 



h(^h(^ lo 



a' 



(3.42) 



The normalized ground state can be expressed in the qubit space as 
-1 



W = 



VW+ 



a' 



{|1100) + |0011) + a|1001) - a|0110)} (3.43) 



The four-partite entanglement in any state \i/j) can be calculated by 
using Eqs. ((lOHl). ( KM} . ( 151251) ) as. 



^. 



\r\ 



1 



where 






/ . ^ijkl 



\ hjik,l^l 



with 



The ground state entanglement can be then calculated to be 

3 ' 



^' ~ (1 + «2: 



l + -a^^a^-l 



(3.44) 



(3.45) 



(3.46) 



(3.47) 
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We plot the four-partite entanglement as a function of U and t (Fig. (2.1a)) 
and as a function of a (Fig. (2.1b)). The entanglement is seen to mono- 
tonically increase as a function of a, saturating at large values of a to the 
maximum value 2. The saturation to the maximum value can be obtained 
either for very large values of U or very small values of t. We can interpret 
this result in the following way: since the total particle number is fixed 
to be 2, the four mode entanglement essentially measures the correlations 
between the spins. The entanglement increases as a function of a because 
the spin correlations increase with a. 




Fig. (3.1): Four-partite entanglement for the Hubbard dimer (at half filling) as a 

)th in 

«(x) = X + -\/l + x^. (a > 1) 



function of (a) U and t (both in energy units) and (b) as a function of a;(-^ ) where 



^4t 



We can also calculate the bipartite entanglement between sites A and B 
using the geometric measure Eqs. ( (15.171 ). (15.181 ). (15.191 )) considering the 
partitions to be {{A t Ai};{Bt B |}} 



K. 



\T\ 



\T\ 



\sep 



where 



and 1 1 T 1 1 sep — 



\T\ 



d{d-l] 



15 /(i\^ 



(3.48) 



(3.49) 



m/2 



. Here As are the generators of SU{4:), there are 
m = 2 partitions (same as the number of sites) and each partition has 
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dimension d 



4. This leads to an intersite entanglement of the form 
2 



Es = 



Vl3a* + 34^2 + 13 - 6 



,, V ^^^ , ^.^ , .^ ^ (3.50) 

The bi-partite entanglement between sites A and B was calculated earlier 
using the von-Neumann entropy 



E 



1 



VN 



log, [2{1 + a^)] -aHog, 



a 



2(l + a2^ 



(3.51) 



(1 + a^] 

We plot the intersite entanglement (the von-Neumann entropy is also plot- 
ted for comparison) as a function of a in Fig. (3.2). It is seen that both 
measures show qualitatively similar behavior, i.e. a monotonically decreas- 
ing entanglement as a function of a saturating at very large values of a. 
The entanglement between the sites A and B decreases as a function of 
a because with increasing on-site repulsion U, the four dimensional lo- 
cal state space at each site gets reduced to a two dimensional local state 
space [EB] due to a suppression of charge fluctuations or in other words, as 
a — > oo the SU{4:) SU{4:) partition goes over to a SU{2) SU{2) par- 
tition. We have explicitly checked that the entanglement obtained in the 
a — > oo limit matches with that obtained for the SU{2) <S) SU{2) partition. 




10 12 



a 




Fig. (3.2) : (a) The bi-partite entanglement between sites A and B calculated with the 
geometric measure as a function of a (a > 1) for the Hubbard dimer at half fiUing. (b) 
The corresponding von-Neumann entropy as a function of a (a > 1). 

We can also discuss bi-partite entanglement with unequal partitioning. 
Consider four modes partitioned into two subsets containing one and three 
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modes respectively. This would correspond for example to one observer 
controlling a register which measures the occupancy of the spin up at site A 
and the other observer controlling a register which measures the occupancy 
of the spin down at site A as well as the occupancy at site B. The bi-partite 
entanglement in this can be obtained as discussed in Sec. (13. 3.21 ). In the 
present case, the two partitions have dimension as di = 2 and d2 = 8 
respectively. This gives, for the entanglement. 



E 



A 



3 63 



j=i k=i 



2 
jk 



28 



(3.52) 



or, using \iIj) in Eq. (13.431 ). 



E = 1.6367 



(3.53) 



Interestingly, the entanglement is independent of a. It turns out that this is 
the maximum value possible for the entanglement (Eq. (I3.52l ))( we checked 
this by maximizing the entanglement given by Eq. (13.521 ) as a function of 
the coefficients in the general state Eq. (13.321 ). We have also checked that 
von-Neumann entropy in this case is also independent of a and has the 
maximum possible value, i.e. 2. 



3.5.2 Three electrons on three sites. 

We next consider the Hubbard trimer, i.e. electrons on three sites with the 
sites A,B,C with periodic boundary conditions. Ignoring the chemical 
potential , the Hamiltonian is 



H(I3) = -t 



E 

j=A,B,C 



^ t it 



+ h.c.) - (3 ^ hj^hj:^ 



j^A,B,C 



(3.54) 



u 



where t > is the hopping parameter and (5 = j. We begin our analysis as 
earlier by mapping from the fermionic to the six qubit space corresponding 
to (^4 t A ], B ^ B I C t C \). In the qubit space, the basis respecting 
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number super-selection rule is given by 

{\k)} k = 7, 11, 13, 14, 19, 21, 22, 25, 26, 28, 35, 37, 38, 41, 42, 44, 49, 50, 52, 56 

(3.55) 
where the twenty basis states are labeled by the six bit binary represen- 
tation of the k values. We numerically diagonalize the Hamiltonian H{f3) 
in Eq. ( 13.541 ) in the basis given by Eq. (13.551 ) at different values of (3 as (3 
increases from zero. The total spin quantum number S and the z compo- 
nent of the total spin Sz commute with the Hamiltonian and can therefore 
be used as good quantum numbers to characterize the states. The ground 
state has total S value = 1/2. The triangular geometry of the three site 
model (with periodic boundary conditions) also leads to an additional sym- 
metry under reflection about one of the medians of the triangle. This leads 
to a two fold degeneracy for the ground state (for a flxed Sz value). Since 
these symmetries are preserved even in the presence of the interaction U, 
the ground state remains two-fold degenerate for all /3 values. 

The entanglement in any state can be calculated in a similar manner as 
shown previously (Sec J3.5.1l ). We now show the results of our calculations 
for the entanglement in one of the ground states. 

The six-mode entanglement as a function of the interaction parameter 
P is shown in Fig (3.3a) while the tripartite entanglement (between the 
sites A,B,C) is shown in Fig. (3.3b) as a function of /3. The tripartite 
entanglement between the sites A^B^C is seen to decrease with increasing 
13 - we interpret this result in a similar way as that for the dimer as due 
to the fact that the local state space at each site decreases with increasing 
f3. We flnd that the six-mode entanglement increases with f3, saturating at 
large values, however, the behavior is not monotonic. Such non- monotonic 
behavior as a function of (3 is also shown by the bipartite entanglement 
between sites A and BC calculated using the geometric measure as well as 
the von-Neumann entropy (Figs. (3.3c) and (3.3d)). 

We can understand the non-monotonic behavior in the following way: 
the bi-partite entanglement measures the entanglement between the the 
sites A and the sites BC. With small increase in (3 from zero, there is an 
increase in the spin fluctuations between B and C which shows up as an 
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initial increase in the entanglement between A and BC. However, with 
further increase in /3, the charge fluctuations get completely suppressed 
leading to an asymptotic behavior similar to that for the dimer. One can 
understand the non-monotonic behavior of the six-mode entanglement in 
a similar way. For small jS, there is a larger correlation between two sites 
which leads to a decrease in the overall entanglement, however with further 
increase in 13, the contribution to the entanglement is solely due to spin 
fluctuations which increases with /3 leading to the observed increase in the 
entanglement as well. 
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Fig. (3.3): (a) Six-partite entanglement for the Hubbard trimer as a function of /3 at 
half filling, (b) Tripartite (site) entanglement as a function of /3 for the Hubbard trimer 
at half filling, (c) The bipartite entanglement between site A and sites BC (see text) as 
a function of /3 for the Hubbard trimer at half filling and (d) The corresponding von- 
Neumann entropy as a function of /3 for the Hubbard trimer at half filling. 



We have also calculated the upper bounds for the entanglements in vari- 
ous partitions in the same way as in Sec J3.4[ These are E"™"^ = 4.42218, 



^max ^ 4.15105, ^™«^ = 6.08767. 
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3.6 Summary and comments 

We have proposed a multipartite entanglement measure for A^ fermions 
distributed over 2L modes (single particle states). The measure is defined 
on the 2L qubit space isomorphic to the Fock space for 2L single parti- 
cle states. The entanglement measure is defined for a given partition of 
2L modes containing m > 2 subsets, using the Euclidean norm of the m- 
partite correlation tensor in the Bloch representation of the corresponding 
multi-mode state, viewed as a m-partite state (see sec. IV and V). The 
Hilbert spaces associated with these subsets may have different dimen- 
sions. The quantum operations confined to a subset of a given partition 
are local operations. This way of defining entanglement and local opera- 
tions gives us the flexibility to deal with entanglement and its dynamics 
in various physical situations governed by different Hamiltonians. We note 
that the concept of locality for indistinguishable fermions is distinct from 
that for distinguishable particles. In the latter case, locality applies to 
spatially separated subsystems. However, spatially separated fermions be- 
come distinguishable. We have shown, using a representative case, that 
the geometric measure is invariant under local unitaries corresponding to a 
given partition. As an application, we have also considered some correlated 
electron systems and demonstrated the use of the multipartite measure in 
these systems. In particular, we have calculated the multipartite entan- 
glement in the Hubbard dimer and trimer (at half filling). We find that 
the bipartite entanglement between sites as computed with the geometric 
measure has a qualitatively similar behavior as a function of the inter- 
action as that of the conventional von-Neumann entropy. We have also 
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calculated the four(six)- partite and the two (three) site entanglement for 
the Hubbard dimer(trimer). We find that the multi-partite entanglement 
gives complementary information to that of the site entanglement in both 
the cases. 

Although the entanglement measures given in this chapter have been 
mainly applied to the study of the multipartite entanglement structure in 
Hubbard dimers and trimers, these measures are completely general and 
can be applied to other fermionic systems. We have also shown jSSl llOOj 
that, viewed as a measure on qubit space, this measure has most of the 
properties required of a good entanglement measure, including monotonic- 
ity. To the best of our knowledge, this is the first measure of multipartite 
entanglement in fermionic systems going beyond the bipartite and even the 
tripartite case. Further, in this chapter we have restricted to applications 
involving ground states of 2L mode systems with 2L < 6. It is straightfor- 
ward to extend the calculations for large number of modes (2L > 6) except 
for the length of the computation. We have shown earlier that jMl llQOj 
for antisymmetric states, the computational complexity of computing the 
measure goes polynomially with the number of parts of the system (here, 
the number of partitions of 2L modes). It would be interesting to extend 
these calculations to larger system sizes where new and interesting results 
might be expected. 
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Chapter 4 

Thermal quantum and classical 
correlations in two qubit XX model 
in a nonuniform external magnetic 
field 



3TT5TT STT3T<HcJ«kiJ|U|l 




A Sanskrit saying (Subhashit) 

People are bound by an amazing chain called hope. Those bound by it keep running. Others stay put as if 

disabled. 



4.1 Introduction 

We have seen in chapter I that composite quantum systems can be in a 
class of states, caUed entangled states, in which the correlations between 
the constituents of the system cannot be achieved in a classical world |120l 



Although all pure entangled states possess such nonlocal quantum 
correlations, there are mixed entangled states which do not, in the sense 
of violating Bell inequalities [3]. The entanglement in quantum states 
and the resulting nonlocal quantum correlations form an area of intense 
research, due to their huge technological promise, especially in the areas 
of quantum communication and cryptography jl24j . However, quantum 
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correlations breaking Bell inequalities need not account for all quantum 
correlations in a composite quantum system in a given state. In order to 
account for the quantum correlation in a given state, we must find some 
means to divide the total correlation into a classical part and a purely 
quantum part. This is particularly important for mixed states, since their 
quantum correlations are many a time hidden by their classical correlations 
{CC). An answer to this requirement is given by quantum discord {QD)^ 
a measure of the quantumness of correlations introduced in Ref. fTB] . 
Quantum discord is built on the fact that two classically equivalent ways of 
defining the mutual information turn out to be inequivalent in the quantum 
domain. In addition to its conceptual role, some recent results [7^, suggest 
that quantum discord and not entanglement may be responsible for the 
efficiency of a mixed state based quantum computer. We believe that QD 
will turn out to be a very useful tool to analyze mixed state quantum 
correlations and their consequences, as mixed state entanglement is very 
difficult and eluding to deal with jl06j . To realize this hope we need a 



viable relation between mixed state entanglement and quantum discord 
[Bl]. The pointers towards such a relation may be obtained by studying 
these properties for various quantum systems. 

Motivated by these considerations, we present , in this chapter, the re- 
sults of our investigation of the amount of QD and CC in a two qubit 
Heisenberg XX chain at finite temperature subjected to nonuniform exter- 
nal magnetic fields Bi and B2 acting separately on each qubit. We study 
two distinct cases namely Bi = —B2 (nonuniform field) and Bi = B2 
(uniform field). In each case, we obtain the dependence of QD,CC and 
entanglement (EOF) in the system on the external magnetic field and tem- 
perature. Such a model is realized, for example, by a pair of qubits (spin 
1/2) within a solid at finite temperature experiencing a spatially vary- 
ing magnetic field. Such Heisenberg models can describe fairly well the 
magnetic properties of real solids |125j and are well adapted to the study 



of the interplay of disorder and entanglement as well as of entanglement 
and quantum phase transitions |126l , 1127] . The variation of entanglement 



1128] and QD J1291 IT^Dj in a two qubit Heisnberg XX chain with external 



magnetic field is already reported. 



4.2 The thermalized Heisenberg system. re 

In order to quantify entanglement in a thermally mixed two qubit state, 
we use concurrence given by [ jl31j 



C = max{Ai - A2 - A3 - A4, 0} (4.1) 

where \i{i = 1,2, 3, 4) are the square roots of the eigenvalues of the oper- 
ator pp in descending order 

p={a\®al)p\al®al), (4.2) 

with Ai > A2 > A3 > A4, and p is the density matrix of the pair qubits; a\ 
and a 2 are the normal Pauli operators. The entanglement is related to the 
concurrence by 

where h{x) = — a:log2X — {1 — x) log2(l — x). Henceforth, in this chapter, 
we denote the entanglement of formation [EOF) by EN. The concurrence 
C = corresponds to an unentangled state and C = 1 corresponds to a 
maximally entangled state. 



4.2 The thermalized Heisenberg system. 

The model Hamiltonian we study is given by 

H = J{S^S^ + ^f^l) + BiS'i + B2SI, (4.3) 

where S" = cr"/2, [a = x, y, z) are the spin 1/2 operators, a" are the 
Pauli operators and J is the strength of Heisenberg interaction. Bi and B2 
are external magnetic fields. As stated in the introduction, by changing 
Bi and B2 separately, we want to study the effects of magnetic field on 
the thermal QD, CC and EN in a very general way. The eigenvalues and 
eigenvectors of H are 

H\m) = -{Bi + B2)m, 
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//|ll) = (5i + 52)|ll), 

H\ij^) = ±D|^±), (4.4) 

where D^ = {Bi - B^f + J^ and |^±) = ]^[|01) + ^^^'^^^ llO)]. We 
denote the eigenvalues corresponding to |00), |11), \il)^) by £"00, £^11, E± re- 
spectively. In the standard basis, {|00), |01), |10), |11)}, the density matrix 
p{T) is given by 

~ ui 

Wi V 

W2 



piT) = i 



where 



U\ = e 



{Bi+B2)/kT 



U2 = e-(^i+^^)/^^, 



wi = cosh( 
W2 = cosh( 






f 



H 7; smhf— -), 

(Bi - B,) D 
7^ smhi-— -), 

_Jsinh(^) 
D ' 



(4.5) 



(4.6) 



and Z = rr[exp(-^)] is the partition function. In the following we select 
I J| as the energy unit and set /c = 1. 



4.3 Quantum Discord. 



We suummarize here the concept of quantum discord, which already dis- 
cussed in chapter 1. In classical information theory (CIT) the total correla- 
tion between two systems (two sets of random variables) A and B described 
by a joint distribution probability p{A^ B) is given by the mutual informa- 
tion (MI), 

I{A, B) = H{A) + H{B) - H{A, B), (4.7) 
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with the Shannon entropy H{p) = — ^jPjlog2Pj- Here pj represents the 
probabihty of an event j associated to systems A,B, 01 to the joint system 
AB. Using Bayes's rule we may write MI as 

I{A,B) = H{A)-H{A\B), (4.8) 

where H{A\B) is the classical conditional entropy. In CIT these two ex- 
pressions are equivalent but in the quantum domain this is no longer true 
|76lll32j . The first quantum extension of MI, denoted by /(/?), is obtained 
directly replacing the Shannon entropy in Eq. (14.71 ) with the von Neumann 
entropy, S{p) = —Tr{plog2p), with p, a density matrix, replacing probabil- 
ity distributions. To obtain a quantum version of Eq. (14.81 ) it is necessary 
to generalize the classical conditional entropy. This is done recognizing 
H{A\B) as a measure of our ignorance about system A after we make a 
set of measurements on B. When 5 is a quantum system the choice of 
measurements determines the amount of information we can extract from 
it. We restrict ourselves to von Neumann measurements on B described 
by a complete set of orthogonal projectors, 11^, corresponding to outcomes 

J- 

After a measurement, the quantum state p changes to pj = [{I^Ilj)p{I^ 
Ilj)]/pj, with / the identity operator for system A and pj = Tr[{I^Ilj)p{I^ 
lij)]. Thus, one defines the quantum analog of the conditional entropy as 
*S'(p|{nj}) = ^jPjS{pj) and, consequently, the second quantum extension 
of the classical MI as [76] 

J{p\{Ilj}) = S'(/) - S{p\{Il^]). The value of J{p\{Il^]) depends on 
the choice of {IIj}. 

Henderson and Vedral [76] have shown that the maximum of j7'(/3|{Hj}) 
with respect to {Hj} can be interpreted as a measure of classical correla- 
tions. Therefore, the difference between the total correlations I{p) and the 
classical correlations Q{p) = supiji^J{p\{Iij}) is defined as 

D{p) = I{p)-Q{p), (4.9) 



giving, finally, a measure of quantum correlations jTSl 1129] called quantum 



discord {QD). For pure states QD reduces to entropy of entanglement 
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highlighting that in this case all correlations come from entanglement. 
However, it is possible to find separable (not-entangled) mixed states with 
nonzero QD [ I76[ 1132] , meaning that entanglement does not cause all non- 



classical correlations contained in a composite quantum system. Also, QD 
can be operationally seen as the difference of work that can be extracted 
from a heat bath using a bipartite system acting either globally or only 
locally 
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4.4.1 Case I: Bi = -B2, and (J > 0). 



In this case, \il)^) is the ground state with eigenvalue E^ = —y/4Bf + J^. 
Other eigenvalues are {0,0, y^'^Bf + J^} for eigenvectors {|00), |11), |^^)} 
respectively. In this case the variation of QD, CC and EN with Bi at 
different temperatures (T=0.2, 0.9, 1.5) is depicted in Figs. (4.1), (4.2) 
and (4.3), respectively. We observe that QD and CC in the thermal state 
coincide for all values of Bi for different temperatures (T=0.2, 0.9, 1.5). 

In order to understand this observation, we take a close look at the 
thermal state. At temperature T, the thermal state is given by 

p = i[|00)(00| + |11)(11| + eV^i^^/T 

|^-)(V^-| + e-V45?+JVr|^+^^^+|]_ (4_;l0) 

This p has the Bloch representation [99] 

1 ^ 

/3= -[/(g)/+^Q(J,(g)(J,], (4.11) 

where Oi (? = 1, 2, 3) are the one-qubit Pauli operators. 

The class of mixed states as in Eq. (14.111 ) have equal classical and quantum 

correlations (like bipartite pure states [133] ) provided 



Q Cj ^ Cfc 
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and 



Ck = -c] (4.12) 

where i ^ j ^ k E {1,2,3} (see the Appendix to this chapter). Here 
Ci^Cj^Ck are the diagonal elements of the correlation matrix defined by 
Cij = Tr{p(Ji <S> (Jj). It is straightforward to check that the thermal state p 
in Eq. ( 14.101 ) which has form of Eq. (14.111 ) satisfies conditions (12). This 
explains the observations in Figs. (4.1), (4.2) and (4.3), that the two qubit 
thermal state p for Bi = —B2 in Eq. ( 14.101 ) gives rise to equal QD and 
CC for all values of Bi and temperature. In order to see why the common 
curve for QD and CC peaks at 5i = we can maximize the expression 
for QD = CC with respect to Bi and check that the maximum occurs at 
Bi = 0. 



0.8 
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0.4 



0.2 



Fig. (4.1) : QD and CC (dashed line) and EN (solid line) as a function of external 

magnetic field Bi = —B2 at T = 0.2 




From Figs. (4.1), (4.2) and (4.3), we also see that EN as a function 
of Bi has a peak at 5i = for T = 0.2, has a dip for T = 0.9 and goes 
to zero over an interval symmetric about Bi = for T = 1.5 [ jl28j . From 
Eq. ( 14.101 ) we see that the concurrence of the thermal state is governed by 
the admixture of the |?/^+) and \il)^) states. We expect concurrence to fall 
as the state \il^'^) classically mixes more and more with the ground state 
\il^^)- For fixed J = 1 and a fixed temperature T, this happens for Bi = 0. 
That is why we got a dip in the EN curve at Bi = 0. The size of this dip 
increases with temperature. In fact the dip touches the Bi axis when, at 
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Fig. (4.2) : QD and CC (dashed line) and EN (solid line) as a function of external 

magnetic field Bi = —B2 at T = 0.9 




Fig. (4.3) : QD and CC (dashed line) and EN (solid line) as a function of external 

magnetic field Bi = —B2 at T = 1.5 



Bi = the concurrence is zero. To see this we note that the concurrence 
for p in Eq. (14.101 ) is given by |134j 



C = — maxjl-yl 



^UiU2, 0} 



(4.13) 



where v,ui,U2 are given in Eq. (14.61 ). Therefore, for Bi = and J = 1, 
C > provided sinh ^ > 1 or T < 1.1346. For T = 1.5 and J = 1, using 
the requirement sinh jt = D, we can find the range of Bi around 5i = in 
which C = 0. This is -1.1456 < Bi < 1.1456. Figs. (4.1), (4.2) and (4.3) 
confirm the corresponding behavior of EN. 

Fig. (4.4) shows the variation of EN, QD, CC with temperature at fixed 
values of Bi = —B2. As expected we have QD = CC for all temperatures. 
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Both EN and QD (CC) curves have plateau at low temperatures corre- 
sponding to their ground state values, as at these temperatures, the ground 
state is not thermally connected to other exited states. Other interesting 
observation is the vanishing of concurrence at a finite critical temperature 
Tc which increases with Bi value, while QD and CC asymptotically go to 
zero with temperature. The increase in T^ with Bi jl28j can be understood 
from the thermal state Eq. ( 14.101 ) which says that higher temperatures are 
required to get a given admixture of \i{j^) and \i/j~^) for higher Bi values. 
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Fig. (4.4) : QD and CC (dashed line) and EN (solid line) as a function of the absolute 
temperature T for (a) Bi = —B2 = 1 (b) Bi = —B2 = 2 



4.4.2 Case II: 5, 



-aBi and (J > 0). 



We now deal with the case B2 = —aBi, a^\ and positive. 5i = satisfies 
both, B2 = —aBi and B^ = — i^i, so that QD = CC at B^ = 0, for all 
temperatures T. For a fixed temperature T, it turns out that QD > CC for 
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Bi ^ if a > 1 and CC > QD for 5i 7^ if < a < 1. This is depicted in 
Fig. (4.5), for a = 2 and a = 1/2 for T = 1.5. The dominance of QD over 
CC (or vice versa) varies continuously with a. This observation gives us the 
key to control the contributions of QD and CC to a two qubit thermal state 
in Heisenberg model via the continuous variation of the applied magnetic 
field. The behavior of concurrence in this case can be analyzed in a way 
similar to the case B2 = —Bi, (a = 1). From Fig. (4.5), we see that for the 
same temperature, the range over which concurrence vanishes depends on 
a, this range decreases monotonically with a. Also, the peak position of 
concurrence (or, EN) on the Bi axis shifts monotonically towards Bi = 
as a increases. Thus the main entanglement features can be controlled by 
varying external magnetic fields. 




Fig. (4.5) : QD (dashed line), CC (dash-dotted hne) and EN (sohd hne) as a function 
of external magnetic field Bi ai T = 1.5 (a) B2 = —2Bi (b) B2 = —Bi/2. 



4.4 Results and discussion. 



83 



4.4.3 Case III: Bi = B2. and (J > 0). 

For uniform external magnetic field Bi = B2, Figs. (4.6), (4.7) and 
(4.8), show the variation of QD,CC and EN with Bi for temperatures 
T = 0.2, 0.9, 1.5, respectively. We see that all three quantities are symmet- 
ric about Bi = where they have their maximum. Further, QD > CC, 
except at Bi = 0, where QD = CC. For higher temperatures, the quali- 
tative behavior of QD and CC remains the same, while EN curve drops 
down below those of QD and CC. This can be qualitatively understood by 
looking at the thermal state given by 

;9=i[e2^^/^|00)(00| + e-2^^/^|ll)(ll| 

+ e^/^|^-)(^-| + e-^/^|^+)(^+|]. (4.14) 

For small temperatures, the entanglement of the thermal state is largely 
dictated by that of l?/'^) and becomes dominant. At higher temperatures, 
admixture due to other states reduces the entanglement, so that QD and 
CC dominate. Such a complementary behavior of entanglement and dis- 
cord can serve as a pointer to wards a possible connection between them. 




Fig. (4.6) : QD (dashed line), CC (dash-dotted line) and EN (solid line) as a function 
of external magnetic field Bx where Bi = B2 aX T = 0.2 



Fig. (4.9), shows the variation of QD,CC and EN with temperature 
for Bi = B2 = I and Bi = B2 = 2. We see that for high temperatures, 
QD hugely dominates EN, showing the robustness of QD with tempera- 
ture. As temperature becomes large QD and CC converge to wards each 
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Fig. (4.7): QD (dashed line), CC (dash-dotted hne) and EN (sohd hne) as a function 
of external magnetic field Bi where Bi = B2 aX T = 0.9 




Fig. (4.8) : QD (dashed line), CC (dash-dotted line) and EN (solid line) as a function 
of external magnetic field Bi where 5i = i?2 at T = 1.5 



other. For larger values of Bi this happens at higher temperatures. As 
the temperature increases, all the coefficients in the thermal mixture Eq. 
(14.141) tend to be equal and the thermal state approaches random mixture. 
Thus it seems that QD and CC approach each other as an arbitrary ther- 
mal state approaches a random mixture. Obviously, for random mixture 
p = \{I®I), QD = CC = 0. A quantitative analysis of the relative behav- 
iors of QD and CC with temperature will be very interesting, but possibly 
have to wait for further developments in the theory. 
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Fig. (4.9) : QD (dashed line), CC (dash-dotted line) and EN (solid line) as a function 
of the absolute temperature T for (a) Bi = B2 = 1 (b) Bi = B2 = 2 
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Fig. (4.10) : ENab (solid hue), QD^ (dashed line), ENae (dot line) and QD^ 
(dash-dotted line), as a function of external magnetic field Bi = —B2 at T = 0.9 
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4.5 Monogamy relations between System and envi- 
ronment. 

It will be interesting to connect our results in previous section with the 
monogamy relations between the EN and the classical correlation jl35j of 
two subsystems (qubits) and the environment 

ENab + CC^ = Sa, 

ENae + CC^ = Sa, (4.15) 

and the relation between EN and QD |136j , 



ENab + ENae = QD^ + QD^, (4.16) 

showing us that EN and QD always exist is pairs. Here A,B label the 
qubits and E stands for the environment. We assume that environment 
(heat bath) comprises the universe minus the qubits A and B so that the 
state pabe is a pure state. Since the variation of all the quantities pertain- 
ing to the system AB with Bi and T are obtained form the XX model, 
we can use Eqs. (( 14.151 ). ( 14.161 )) to find the corresponding dependence of 
ENae and QD^ on Bi and T. Figs. (4.10) and (4.11), (for Bi = -B2) 
show the variation of ENae and QD^ with Bi and T. 

The monogamic relations also help us establish a necessary and sufficient 
condition for QD^ = CC^ when the environment is present. This is : 
QD^ = CC^ if and only if \Iab = ENae + ENab - QD^. To prove 
the necessity we note that when QD^ = CC^, (that is, QD^ = ^Iab), 
Eq. ( 14.161 ) can be written as 



1 
2 
Now suppose Eq. ( 14.171 ) is true. Then using Eq. ( 14.161 ) we have 



Iab = ENae + ENab - QD^. (4.17) 



^^^ = 2^A^ 
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which imphes QD^ = CC^. Figs. (4.12) and (4.13), show the variation 
of both sides of Eq. (14.171 ) with Bi and T which estabhshes Eq. (14.171 ) for 
the XX model. 
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Fig. (4.11) : ENab (solid line), QD^ (dashed line), ENae (dot line) and QD^ 
(dash-dotted line), as a function of the absolute temperature T ior Bi = —B2 = 1. 
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Fig. (4.12) : ^Iab (solid line), ENae + ENab — Q^^ (dashed line), as a function of 
external magnetic field B2 = —Bi at T = 1.5 
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Fig. (4.13) : jIab (solid line), ENae + ENab — QD^ (dashed line), as a function of 

the temperature T for B2 = — -Bi = — 1. 



4.6 Summary and comments 

In this chapter, we have studied the variation of. QD, CC and EN in two 
qubit XX Heisenberg chain as functions of independently varied magnetic 
fields Bi and B2 on each qubit and also with temperature. We deal with two 
cases Bi = —B2 (nonuniform field) and Bi = B2 (uniform field). Our first 
observation is the complementary behavior of entanglement and QD/CC. 
For the nonuniform magnetic field, we get the interesting observation that 
QD and CC are equal for all Bi = —B2 values as well as all temperatures. 
Surprisingly, this observation is explained quite simply, using the symmet- 
ric form of the thermal state. A very interesting observation is that the 
relative contributions of QD and CC can be tunably controlled by varying 
the applied magnetic field. Another interesting finding is that the equality 
of QD and CC of the subsystem (qubits) imposes a constraint on the dis- 
tribution of. QD and EN over the subsystem and its environment. Further 
investigation of general Heisenberg models like XXZ along these lines may 
turn out to be interesting and fruitful. 



Appendix: 

Theorem: If the quantum state has the Bloch representation [99] 



1 ^ 

/?=-[/ (g)/ + ^ Cia^ (g) ai\ , 



i=l 



4.6 Summary and comments 



and Q = Cj > Ck and Ck = — cf where i ^ j ^ k E {1,2,3}, then 
this state contains the same amount of quantum and classical correlation 
{QD = CC). 

Proof: In Ref. [i86j S. Luo evaluated analytically the quantum discord 
for a large family of two-qubit states, which have the maximally mixed 
marginal and their Bloch representation is 



/?=-[/ (g)/+^ c^a^ (g) ai\ . 
i=l 

For this class of quantum states the quantum mutual information is 

I{p) = i[(l-Ci-C2-C3)/0^2(l-Cl-C2-C3) + (l-Ci+C2+C3)/0^2(l-Cl + C2 
C3) + (l + Ci-C2 + C3)/05r2(l + Ci-C2 + C3) + (l+Ci + C2-C3)/05r2(l + Cl + C2-C3 

We substitute the conditions above in the quantum mutual information. 
Puting c = ci = C2 > C3 and C3 = — c^, we get, 

I{p) = i[(l-2c+c2)/0(/2(l-2c + c2) + (l-c2)/0(72(l-c2) + (l-c2)/0^2(l- 

c^) + (1 + 2c + c^)log2{l + 2c + c^)] After some algebraic simplification, we 

get 

I{p) = (1 - c)log2{l - c) + (1 + c)/o^2(l + c) 

which equals 2CC as in Ref. [HS]- It is also easy to check that the above 

argument goes through when c = ci = C3 > C2 and C2 = — c^ and when 

c = C3 = C2 > ci and ci = — c^, to get T{p) = 2CC. Thus, 

QD{p)=I{p)-CC{p) = CC{p). 
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Chapter 5 

Geometric measure of quantum 
discord for an arbitrary state of a 
bipartite quantum system. 






Author unknown 

A STATEMENT WHICH HAS AN ESSENTIAL MESSAGE, WHICH IS SHORT, SIMPLE AND UNAMBIGUOUS IN EXPRESSION WHILE 
GENERIC AND UNIVERSAL IN APPLICATION, IS CALLED A PRINCIPLE OR A FORMULA. 



5.1 Introduction 

Quantum discord, (OUiver and Zurek [TB]), is a measure of the discrep- 
ancy between quantum versions of two classically equivalent expressions 
for mutual information and is found to be useful in quantification and ap- 
plication of quantum correlations in mixed states. It is viewed as a key 
resource present in certain quantum communication tasks and quantum 
computational models without containing much entanglement. A step to- 
ward the quantification of quantum discord in a quantum state was by 
Dakic, Vedral, and Brukner [53] who introduced a geometric measure of 
quantum discord and derived an explicit formula for any two-qubit state. 
Recently, Luo and Fu jl37j introduced a generic form of the geometric 



measure of quantum discord. Using these results we find an exact formula 
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for the geometric measure of quantum discord for an arbitrary state of a 
mxn bipartite quantum system. This geometric measure is experimentally 
accessible. 

Correlations in quantum states, with fundamental applications and im- 
plications for quantum information processing, are usually studied in the 
entanglement- versus-separability framework [ j251 1138] . However, entangle- 
ment, while widely regarded as nonlocal quantum correlations, is not the 
only kind of correlation. An alternative classification for correlations based 
on quantum measurements has arisen in recent years and also plays an im- 
portant role in quantum information theory [ J1391 I86l 11401 1141] . This is 



the quantum-versus-classical paradigm for correlations. In particular, the 
quantum discord as a measure of quantum correlations, initially introduced 
by Ollivier and Zurek [76] and by Henderson and Vedral [79] , is attracting 
increasing interest [143-164]. Recall that the quantum discord of a bipar- 
tite state /9 on a system H°^ H^ with marginals p"' and p^ can be expressed 
as 

Q{p) = rmn{I{p) - I{U%p))}. (5.1) 

11" 

Here the minimum is over von Neumann measurements (one dimensional 

orthogonal projectors summing to the identity) 11" = {H^} on subsystem 

a, and 

H» = ^(H^®/V(n^®/^) 

k 

is the resulting state after the measurement. I{p) = S{p"') + S{p^) — S{p) 
is the quantum mutual information, S{p) = —trplnp is the von Neumann 
entropy, and I^ is the identity operator on H^. Intuitively, quantum discord 
may thus be interpreted as the minimal loss of correlations (as measured 
by the quantum mutual information) due to measurement. This formula- 
tion of quantum discord is equivalent to the original definition of quantum 
discord by Ollivier and Zurek [76]. Recently, Dakic et al. introduced the 
following geometric measure of quantum discord [183] : 

D(p) =min||/9-x|P, (5.2) 

X 
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where the minimum is over the set of zero-discord states [i.e., Q{x) = 0] 
and the geometric quantity 

11/3 -XIP --trip-xf 

is the square of Hilbert-Schmidt norm of Hermitian operators. A state x 
on H"^ H^ is of zero discord if and only if it is a classical-quantum state 
|77l [7Bj , which can be represented as 

m 

X = '^Pk\k){k\^Pk, (5.3) 

fc=i 

where {pk} is a probability distribution, {\k)} is an arbitrary orthonormal 
basis for H^ and pk is a set of arbitrary states (density operators) on H^. 
Dakic et al [H3] obtained an easily computable exact expression for the 



geometric measure of quantum discord for a two qubit system, which can 
be described as follows. Consider a two-qubit state p expressed in its Bloch 
representation as (see below) 

1 ^ 

p = - (r /^ + ^(x.a, (g) /^ + r (g) y,a,) 

i=l 
3 

{ai} being the Pauli spin matrices. Then its geometric measure of quantum 
discord is given by [53] 

D{p) = \{\\x\\'+\\T\\'-Xma.). (5.4) 

Here x := (a:i,a:2, ^3)^ and y := (2/i,?/2, 2/3)^ are coherent (column) vectors, 
T = (tij) is the correlation matrix, and Xmax is the largest eigenvalue of the 
matrix xx^ + TT*. The norms of vectors and matrices are the Euclidean 
norms, for example, ||x|p := Xli^^ Here and throughout this article, the 
superscript t denotes transpose of vectors or matrices. 

In order to obtain the exact formula for the quantum discord in an 
arbitrary bi-partite state, we set up the following scenario. Consider a 
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bipartite system H'^ H'^ with dim H"^ = m and dim i7^ = n. Let L(H"') 
be the space consisting of all linear operators on H"^. This is a Hilbert 
space with the Hilbert-Schmidt inner product 

{X\Y) = trX'^Y. 

The Hilbert spaces L(H^) and L(H'^ H^) are defined similarly. Let 
{Xi : i = 1,2,- •• ,m^} and {Yj : j = 1,2, ••• , n^} be sets of Hermi- 
tian operators which constitute orthonormal bases for L[H°') and L[H^) 
respectively. Then 

trXiXi' = 5ii', trYjYji = 5jj'. 

{Xi^Yj} constitutes an orthonormal (product) basis for L{H°'®H^) (linear 
operators on i7° ® H^). In particular, any bipartite state p on H"^ ® H^ 
can be expanded as 

p = ^c,^X,®Yj, (5.5) 

with Cij = tr{pXi ^Yj). 

Quite recently, S. Luo and S. Fu introduced the following form of geo- 
metric measure of quantum discord [ jl37j 



D{p) = tr{CC*) - maxtr{ACC*A'), (5.6) 



A 



where C = [cij] (Eg. (15.51) ) is an m^ x n^ matrix and the maximum is taken 
over all m x m^-dimensional isometric (see below) matrices A = [aki] such 
that 

aki = tr{\k){k\X^) = {k\X,\k), /c = 1, 2, . . . ,m ; z = 1, 2, . . . , m^ 

(5.7) 

and {\k)} is any orthonormal basis in H"^. We can expand the operator 
\k){k\ in the basis {Xi} as 

\k){k\=^ak^X,, /c = l,2,...,m. (5.8) 
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A = [aki] is an isometry in the sense that AA^ = P and the row vectors dk 
of the matrix A satisfy 

11411' = ^4 = 1- (5-9) 

i=l 

Further, using their definition it immediately follows that 

m 

Y^ak^ = trX,. (5.10) 

fc=i 

We need the structure of the Bloch representation of density operators, 
which can be briefly described as follows |121l 1^ 1106] . Bloch representa- 
tion of a density operator acting on the Hilbert space of a d-level quantum 
system C^ is given by 

P = ^ad + I]s,A,), (5.11) 

i 

where the components of the coherent vector s, defined via Eg. (15. Ill ), are 
given by Si = ^tr{pXi). Eg. (15. Ill ) is the expansion of p in the Hilbert- 
Schmidt basis {Id,Xi',i = 1, 2, . . . , d^ — 1} where A^ are the traceless or- 
thogonal hermitian generators of SU{d) satisfying tr{XiXj) = 26ij |121l 



These generators are characterized by structure constants fijk (completely 
antisymmetric tensor) and gijk (completely symmetric tensor) of Lie alge- 
bra su{d) via 

~ ~ 2 

AiAj = -Sijid + ifijkXk + QijkXk (5.12) 

A systematic construction of the generators of SU{d) is known [25] and 
are given by 

{^i}L^^ = {ujk,Vjk,(^i}, (5.13) 

where 

^jk = \j){k\ + \k){jlvjk = -ii\j){k\-\k){j\) l<j<k<d 

^i = Jj^^ili\3){3\-l\l^l){l + l\) l<l<d-l 

(5.14) 
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with {|s)}f^2 being some complete orthonormal basis in Hd. In what fol- 
lows, we take the generators of SU{d) as defined via Eqs. (15. 13115. 141 ). 

We can represent the density operators acting on a bipartite system 
i7° (g) H^ with dim{H°^) = m and dim{H^) = n, as 



P = {Im ® /n + ^ XiXi (g) /„ + ^ yjim ^ Xj 

1 1 LI V 



+ ^^UjK®\j), 



(5.15) 



^j 



where Aj, z = 1, . . . , m^ — 1 and A^, j = 1, . . . , n^ — 1 are the generators 
of SU{m) and SU{n) respectively. Notice that x G M™ '^ and y E W^ ^^ 
are the coherence vectors of the subsystems A and B, so that they can be 
determined locally. These are given by [i96j 



m ~ m ~ 

Xt = —tr{pX^ (g) In) = —tr[pAXi) 

Ti ~ Ti ~ 

Vj = 2^r{plm Xj) = -tr{pBXj), 

where pA = trsip) and pb = tr^ip) are the reduced density matrices. The 
correlation matrix T = [tij] is given by 



T = [Uj] 



mn 



[tr{p\®Xj)]. 



5.2 Main result. 



We find the maximum in Eq. (l5.6p to obtain an exact analytic formula, as 
in the two-qubit case (Eq. (15.41) ) [ 183] . 

Theorem 1. Let /o be a bipartite state defined by Eq. (15.151) . then 



D(p) = 



m?n 



\x\\ -\ — \\l \ 
n 



TO— 1 

1=1 



IV^-I 



(5.16) 



where r/j, j = 1, 2, • • • , m^ — 1 are the eigenvalues of the matrix (xx* + ^^^) 
arranged in increasing order (counting multiplicity). 



5.2 Main result. 
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We prove this theorem for arbitrary (finite) m and n. 
We choose the orthonormal bases {Xi} and {Yj} in Eg. (15 .51 ) as the 
generators of SU{m) and SU{n) respectively fl21j . that is, 



Xi — —^Imj Yl — —^Ifi^ 

'm \/n 



and 



X, = ^ Vi, ^ = 2, 3, . . . , m^ 
Yj = ^^-b J = 2, 3, . . . , nl 



Since trAj = 0; z = 1, 2, • • • , m — 1, we have, via Eg. (15 .101 ). 



m 



y^ flfci = ^^^i = trX^^i = 0, i = 2, . . . , m' 



Therefore, 



771—1 



a^ 



= -^flfcz, ^ = 2,3, •• • ,m' 



(5.17) 



fc=i 



We now proceed to construct the m x m^ matrix A defined via Eg. (15. 71 ). 
We will use Eg. (15. 171 ). The row vectors of A are 



Ofc = {dki, a/c2, • • • , flfcrnO? k = 1,2,. . . ,m. 



Next we define 



m 



Cfc 



V m — 1 
and using Eg. (15. 171 ). we get 



ak2, aki, ..., akm^), k = 1,2, . . . ,m - 1 



(5.18) 



^m 



TO— 1 

k=l 



We can prove 



\ek\\ =1 k = 1,2, ... ,ra 



(5.19) 



(5.20) 
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using the condition ||afc|| = X]j=i a^^ = 1 (Eg. ( 15.91 )) and using Eg. ( 15. 101 ) 
with i = 1, namely, aki = tr{\k){k\Xi) = -j^. Further, isometry of the A 
matrix (74^4* = /) implies 

e,e] = -, jV ^ = 1, 2, . . . , m. (5.21) 

■' m — 1 

We can now construct the row vectors of m x m^ matrix A, using Eg.( l5.18f ) 
and Eg. (Eli). 

1 



ak 



'1, Vm — lefc), /c = 1, 2, • • • , m — 1 



ttm. — 



XW: 



m 



TO— 1 



fc=l 



The matrix A is, in terms of its row vectors defined above, 

/ 1 \/m — lei \ 

1 



A = 



1 



m 



1 



y/m — lei 
\/m — le2 
\/m — les 



\ 1 -Vm-l)^^^! efc y 

We get the elements of C = [cij] = [tr{pXi0Yj)] using the definitions of 
the bases {Xi} and {Yj} given above, in terms of the generators of SU{m) 
and SU{n). These are 

cii = -^, 



" rriy/n " ' ' 



2_ 
TO,n *.? 



where x^ and ?/j are the components of the coherent vectors in Eg. ( 15.151 ) 
and ijj are the elements of the correlation matrix T. The matrix C can 
then be written as 



C = 



'mrm nwm 



V^y^ 



V2 



-X ^-T 

m\/n mn 
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and 






n^rm?- 



■T||l (5.22) 



Having constructed the matrices A and C, we get, for tr{ACC A), 



tT{ACC'A') = -l--r-M?^^^'^ ^ 



y ' 9 



771—1 771—1 777—1 

.j = l 7 = 1 j = l 

(5.23) 



where 



2TT^ 

G = xx* + (5.24) 

n 



is the (m — 1) x (m — 1) real symmetric matrix. For m = n = 2, as in 



triACC'A') = 1 [1 + ||^||2 + eiGe\] , (5.25) 

which attains the maximum of tr{ACC^A^), when e\ is an eigenvector of 
matrix G for largest eigenvalue, which give us Eq. (15^) . The eigenvectors of 
G span R"^ ~^ and form a orthonormal basis of R"^ ^^. Let r/i, 772, ... , r/j„2_x 
be the eigenvalues of G arranged in increasing order (counting multiplic- 
ity). Let (l/i), 1/2), • • • , \fm^-i}) be the corresponding orthonormal eigen- 
vectors of G. Here and in what follows we denote the column vectors in a 
orthonormal basis in M™ ^^ by ket and the corresponding row vectors by 
bra. Whether a ket (bra) vector belongs to M™ ~^ or H°^ can be understood 
with reference to context. 

The last constraint on vectors {cj} making up matrix A is that {cj}, j = 
1, . . . , m must be the coherent vectors of m pure states comprising an or- 
thonormal basis of H°' [HS] • We can span all such sets of coherent vectors 
(each making up an A matrix) using Lemma 1 proved in the Appendix. 
As seen from Lemma 1, in order to span all the sets of coherent vec- 
tors making up matrix A, we can start with the set of coherent vectors 
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{nj},j = l,...,m corresponding to some orthonormal basis in H"^ and 
then obtain all possible other sets of coherent vectors composing matrix 
A via Eq.(Al) corresponding to the unitary transformations U taking the 
initial orthonormal basis to other orthonormal bases in H"^. We choose the 
initial orthonormal basis with coherent vectors {fij} to be the standard 
(computational) basis {|/c)} /c = 1, . . . , m. We can write 



m^ — 1 



hj = ^ e^J\^) J = l,...,m 



(5.26) 



i=i 



where now the set {\i)} i = 1, . . . , rn? — 1 is the orthonormal basis in W^ ~^, 
the ith basis vector consisting of 1 at the ith place and zero at all other 
places, e = [cij] is a m^ x m matrix with its jth column being vector fij 
expressed in {\i)} basis. Obviously, 



^fcj 



m 



2 m 



-^iJl^klJ)' 



Substituting Eq.(Al) and Eq.( 10BI ) in Eq.(E23D we get 



tr{ACC'A') = ^ { - + ^\\y\\^ + 



2{m-l] 



m—\ rri^—l rri^ — l 

E E E ^kjepj{k\OGO'\p) 

j=l k=l p=l 



771—1 771—1 777^ — 1 777,^ — 1 

+ E E E E ^Me,j{k\OGO'\p) 

i=l j=l k=l p=l 



(5.27) 



\m^ — l 



We write G = X^^^ ''lq\fq){fq\^ with its eigenvalues arranged in non- 
decreasing order, giving 



ill 2 2 2(m-l) 



777^ — 1 777^ — 1 m? — l 



tr{ACC'A') = -{- + —\\y\\'^ 



w?n 



m \ n n 

777—1 777—1 777—1 



j-i 



7-1 J-1 



E E E^^(^i^iA)(Ai^'i^) 

fc=l p=l q=l 

(5.28) 
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We substitute tkj = J^^^UMj) in Eq.(|5.28|), to get 



-I \ -I rf -1 m^—l m — 1 m — 1 

tr{ACC'A') = -{- + -\\y\\' + ^ E E E V,{k\0\f,){f\\0'\p) 
m \ n n^ mn ^-^ ^-^ ^-^ 

y fc=l p=\ q=\ 

rn— 1 m—\ to— 1 

-J=l *=1 J=l 



(5.29) 



To evaluate the sums in Eg. ( 15.291 ) we need the average values of Afc 
operators in the standard (computational) basis states in E.^ . Using the 
definition of Afc,A; = l,...,m^ — 1 operators in Eqs. (15. 13115. 141) |95] and 
the fact that \i), |j) in Eg. (15.291) correspond to the standard basis in iif^, 
we see that {i\\h\i) 7^ 0, /c = 1, 2, • • • ,m^ — 1 only for k = (/ + 1)^ — 
1, 1 < / < m - 1. Also, (j|Ap|j) 7^ 0, j9 = 1,2,--- ,m^ - 1 only for 
]7 = (s + 1)^ — 1, 1 < s < m — 1. For A(/+i)2_x or A(s+i)2_i we have, 



/(.) 



'^(/(s) + l)2-l 

where 



/(s)(/(s) + 



Yy(El^')01-^(^)l^(^) + i)(^(^) + i|)- 



j=i 



1 < /(s) < m- 1. 
Using the definition of \(i[s)+\Y-\ above, we can easily prove 



771—1 



E 0' l^a+i)^-i b') (j l^(.s+i)2-i Ij ) 



2 ifl</ = s<m-2 
— if / = s = m— 1 
otherwise 



TO— 1 771—1 

E^^i^(^+i)'-ii^)E^-^'i^( 



,+ l)2_l|j) = i - 



i=l 



2(to-1) -r , -, 

^ 11 / = s = m — 1 

otherwise 







Following the para after Eg. (15.291) and the above eguations. Eg. (15.291 ) can 
be reduced to 
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m^ — 1 



t At^ 



tr{ACC'A 



m \ n n^ mn ^-^ 



TO— 1 



^|((( + 1)^-1|0|/, 



/=1 



(5.30) 

To get the maximum of tr{ACC*A*) over A matrices we use Cauchy- 
Schwarz inequality 

m + If - i|oi/,)p < ((( + ly- -i\(i+ If - i)(/,io'oi/,> = 1, 

which tells us that the desired maximum is obtained by choosing O in 
Eq.( l5.3Q[ ) to be that orthogonal matrix which takes the eigenbasis of G 
matrix to the standard basis in R"^ ~^ We denote the inverse of this or- 
thogonal matrix by 0\^^^. Thus we get, 



maxtrMCC'^O = - J i + ^\\y\\^ + — 
A m \ n n^ mn 



TO— 1 



Yl ^('+1)^ 



1=1 



Finally, Eq.(j02D, Eq.(10B and Eq.(ES!) together imply 



. (5.31) 



D{p) = 



m?n 



\x\ 



n 



TO— 1 



IV^-I 



1=1 



To complete the proof we have to show that the vectors 



^max"^] '■> J ~ -*-' • 



, m 



1 



(5.32) 



are the coherent vectors of some pure states in H"^ given that the vectors 
{nj} are the coherent vectors of the computational basis states in iif". A 
necessary and sufficient condition for the vectors {eH to be the coherent 
vectors of some pure states in H^ is 



J J 3 ■ 



m 



J = 1, • 
where the star product of vectors a and b is defined as 



(5.33) 



^- A /m(m-l) 1 Y^ , 



i,j 
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gijk being the elements of the completely symmetric tensor for SU{m), de- 
fined as in Eq. (l5.12j ). The orthogonal operator O^^^ maximizing tr{ACC^A^] 
is given by 

k 

so that the vectors e^- are, using Eg. ( 15.261 ) 



■m^—\ 



ej = O^^^^^hj = ^ e,j\fi) j = l,...,m. 

i=l 

Thus vectors e* have the same components with respect to the rotated 
basis {|/i)} as the components of fij with respect to the standard basis. 
(That is, both the standard basis and the vectors hj undergo the same 
rigid rotation.) It follows that e* satisfy Eg. ( 15.331 ) if vectors fij do. This 
completes the proof. 

i?emark 1 : 

The choice of {Afc} operators, as defined in Egs. ( 15. 13115. 141 ) [95] is not a 
restriction causing loss of generality. An arbitrary permutation of {Afc} by 
a permutation matrix P replaces G by PGP* = Y17=i '^P{q)\fp{q)){fp{q)\ 
causing appropriate re- arrangement of eigenvectors \fq) and the corre- 
sponding eigenvalues rjq. This leaves Eg. ( 15. 3 If ) invariant under such a per- 
mutation. 

i?emark 2 : 

We expect D{p) to vanish if /? is a classical guantum state as in Eg. ( 15.31 ). 
This follows trivially from the fact that the generators of the SU{m) group 
are traceless, that is, 

m 

tr{\i) = J^i^Mk) = z = 1, 2, . . . , m^ - 1, 

where {|/i;)}; A: = 1, 2, . . . , m is an arbitrary orthonormal basis in H"^. 

We emphasize that the geometric measure of guantum discord in Eg. ( 15.161 ) 
is experimentally accessible, as all the guantities involved depend on the 
average values of {Aj} and Aj0 Aj (generators of SU{m) and SU{n) respec- 
tively) which are hermitian operators and hence are in principle measurable 
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for any m and n. However, it turns out that feasible experimental meth- 
ods to measure A^, Xj are available only for qubits and qutrits fl64l 1165] in 
which case they are simply related to spin angular momentum components 
of a spin-^ and spin 1 particle respectively. Thus we can experimentally 
determine the quantum discord in a bipartite system comprising qutrits 
and qubits. It is important to note that this experimental determination 
of the quantum discord does not require a detailed knowledge of the state 
of the system and can apply to any stage of quantum information process 
in which the state of the quatum system may not be known. In general, 
determination of the unknown state of a quantum system is a formidable 
task. 

For m = n = 3, we have 



Dip) = 



27 



W\\' + ^\\T\\'-{m + m) 



(5.34) 



We now give some examples using Eq. (l5.16f ). 

Example 1. We consider the m x m-dimensional Werner state 

m- z ^ mz - 1 ^ r..i 

P = — — -^ + — — -F. ^ e [-1, 1] 



rw — m 



171'^ — m 



with F = ^f^i \k){l\ \l){k\- First, we calculate the coherent vector x. We 
have 



m 



m 



^i = —tr{p\i®I) = — 

+ 



m — z 
m^ — m 
mz — 1 



tr{\ (g) /) 



rri'^ — m 



tr{F{X,(^I)) 



(5.35) 



Since tr{\i) = 0, (i = 1, 2, . . . , m^ — 1), the first term in the above equation 
is zero. This gives 



Xi 






mz 



m^ 



m ^ — ' 



m 
~2 



mz 



1 



m- 



m 



tr{\ 



0, 



so that X = 0. 
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Now, we calculate the elements of the correlation matrix 



m 



tij = —-tr{p\i (g) A 



m 



'], 



m — z 



m'^ — m 



tr{Xi Xj 



'mz — 1 , ^,~ ~ ,, 

+ —. tr{F{X, ^ Xj)) 

rri'^ — m 

The first term in the above equation is zero, so that 



(5.36) 



m 



mz — 1 



m-' — m 



j2{i\Vxj\i) 



m 



mz — 1 



m'^ — m 



tr{XiXj 



Using the orthonormality of the generators {A^} we get 



m 



nj 



mz — 1 



m'^ — m 



(2^. 



u; 



or, 



l^ii 



m^ (mz — 1) 



2 (m^ — m) 
Thus T matrix is diagonal and all diagonal elements are equal. This gives. 



2iiri 



m^ — 1 



m 



„ ,,1, i ^ 

m < ^ m 



m 



i=l 



m^ (mz — 1)' 
4 (m"^ — m) 

m.{wiz — 1)^ 



m 



m 



1) 



2(m2-l) ' 
the eigenvalues of ^TT* are all equal : 

m,^{m,z — 1)^ 
2(m^ — m)^ 



^z 



Substituting \\x\ 



and the expressions for 



2||r| 



and {iji} as above 



in Eq. (l5.16j ) we get, after some algebra, for the quantum discord of the 
7TT. X ?7T--dimensional Werner state 



D(p)^ 



[m,z — 1)' 



m.{m. — l)(m + 1)^ 
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A Werner state is separable if and only if z G [0, 1], but its geometric 
measure of quantum discord vanishes if and only if. z = 1/m. 

Example 2. We consider the particular case m = n = 3 and the bipartite 
pure state 

|^> = i|ll) + i|22) + i=|33), 

then by straightforward evaluation based on the original definition, its 
geometric measure of quantum discord is D(|?/^)(?/^|) = |. We show that 
Eq. (l5.16l ) yields the same result. 

The detailed calculation is as follows. 



\^){^\ = lU^ 



I 



V3 



4 



/0As 



Vs 



A, 



9 



9 

A1 + --A2 



9 



--Xa X4 



9 



:A5 A5 + 



9 



4^2 4V2 " ^ 4x/2 

Using this expression for \i/j){i/j\, we get 



Afi Xf 



9 



- 9^ 

A2 + 0A3 

o 

15 



A, 



4^2 



A7 A7 + --Ag Ag} 



G = ff* + -TT^ 
3 



/ 27 
' 32 













27 
32 














27 
32 













27 
16 













27 
16 













27 
16 













27 
16 





\ 













32 / 



and from Eq. ( 15.341 ) we get, using the values of ||x| 



12 2| 



Dip) = ^ 



3 2 279 

h - X 

16 3 16 



'27 81^ 
32 ^32 



|r|p and 7/3 + T/g, 
5 



Example 3. We consider the two qutrit state 

p = p\e){e\ + (1 -j?)- 



(5.37) 



where |e) = ^(|2)® |2) + |3)® |3) + |2)(^|1) + |1)(^|2) + |1)(^|3) + |3)® |1)), 



V& 



I is the identity operator and {\i)] i = 1,2, 3} is the standard basis in C^. 
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Fig. (5.1) shows the variation of D{p) (Eg. (15. 161 )) and the lower bound on 
D{p), as given in fl37j . namely, triCC^) — YlT^i^i (where {Aj} are the 
eigenvalues of CC^ listed in the decreasing order, counting multiplicity), 
with p. We see that D[p) (Eg. (15. 16] )) dominates this lower bound for p > 
0.2. 



0.7 



0.5 



0.3 



0.1 ^ 




Fig.(5.1) 



Quantum discord (Eq.f l5.16"|) ) (solid line) and its lower bound (Eq.(5.6) in 
[137] ) (dashed line) as a function of p (Eq. (l5.37p ). 



Example 4- We consider the two gutrit state 

p =;p|ei)(ei| + (1 -p)\e2){e2\ 



(5.38) 



where |ei) = i|ll) + i|22) + ^|33) and |e2) = ^(|2) ® |2) + |3) ® |3) + 
|2) (g) |1) + |1) (g) |2) + |1) (g) |3) + |3) (g) |1)). Fig.(5.2) shows the variation 
of D{p) (Eg. (15. 16] )) and the lower bound on /^(p), as in jl37j . namely, 
triCC^) — YlT=i ^i (where {A^} are the eigenvalues of CC^ listed in the de- 
creasing order, counting multiplicity), with j9. We see that D{p) (Eg. (15. 16] )) 
dominates this lower bound. 
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0.7 



0.5 



0.1 



Fig. (5.2) : Quantum discord (Eg. (15. 16"]) ) (solid line) and its lower bound (Eq.(5.6) in 
[137] ) (dashed line) as a function of p (Eq. (l5.38p ). 
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Appendix to the chapter 5 



A CHILD OF FIVE WOULD UNDERSTAND THIS. SeND SOMEONE TO FETCH A CHILD OF FIVE. 

Groucho Marx 

In this appendix we prove the lemma 1. For a given pure state \i){i\ and 
an unitary operator U acting on i7°, there exists an orthogonal operator 
O = [Oais] acting on MJ^ ~^ such that 

e, = h]0 (Al) 

where hj is the coherence (column) vector of the state \i){i\ and Cj is 
the coherence (row) vector of the state U\i){i\W. 

Proof: We use the following fact which is easily proved. 

UXaU^ = J2 Oap>^p a = l,...,m^-l. {A2) 

P 
for some orthogonal operator O = [Oap] ■ 
We now have, using Eq.(A2) 

^(e,)/3A^ = U\i){i\U'^ = J2inj)o.UKU^ = ^ ( ^(n,).0,J A^ 

/3 a 13 \ a J 

which proves the Lemma 1. 
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Chapter 6 

Summary and Future Directions 






The first chapter is introductory and does not present any new work. 
Therefore, we start summarizing from the second chapter. 

In chapter 2, we deal with the problem of controlled production of entan- 
glement that is, producing the multipartite quantum systems in required 
entangled states. The motivation of this work is as follows. It is now well 
known that entanglement is a resource for accomplishing various kinds of 
quantum information processing tasks and quantum communication proto- 
cols, thus future technology may require uninterrupted supply of entangled 
quantum systems just as the present day technology requires uninterrupted 
supply of energy. Therefore, it is worth searching for the general methods 
of preparing entangled quantum states of multipartite quantum systems. 
This can be achieved, at least in principle, because different interaction 
Hamiltonians can be expressed in a single generic form. We have imple- 
mented some protocols for entanglement generation for two qubits, two 
qutrits and three qubits using the geometric measure of entanglement sug- 
gested by us pifTOD]. 
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A possible question that can be addressed is the effect of ancillas in en- 
tanglement production. In ref. [42j ancillas are proved to improve efficiency 
of entanglement production in the two qubit case. 

A possible analytical approach to find the entanglement capacity for 
two qutrits and three qubits is another challenge. 

In chapter 3, we have given a geometric measure of entanglement for a 
multipartite fermionic system. We construct a isomorphism which enables 
us to view a fermionic (Fock) state as multiqubit state to which we can 
apply the geometric the geometric entanglement measure invented by us 
[ 1991 llOOj . We think that our way of defining the local operations in the case 
of fermionic systems has removed the ambiguities prevailing before. We 
have tested our measure on the generically important systems as Hubbard 
dimer and trimer. 

The basic open question in this area is the possible role of entanglement 
in the wide range of cooperative phenomena in many body systems, in 
particular quantum phase transition. We can use our measure to quanti- 
tatively study the dependence of entanglement on various interaction pa- 
rameters in the model Hamiltonian, which in turn are connected to the 
physical behavior of systems. One of the immediate problems is the fate 
of entanglement when the Coulomb interactions between different sites are 
different (in Hubbard trimer). 

In chapter 4, We investigate how thermal quantum discord {QD) and 
classical correlations {CC) of a two qubit one-dimensional XX Heisen- 
berg chain in thermal equilibrium depend on temperature of the bath as 
well as on nonuniform external magnetic fields applied to two qubits and 
varied separately. We show that the behavior of QD differs in many un- 
expected ways from thermal entanglement {EOF). For the nonuniform 
case, {Bi = —B2) we find that QD and CC are equal for all values of 
{Bi = —B2) and for different temperatures. We show that, in this case, 
the thermal states of the system belong to a class of mixed states and satisfy 
certain conditions under which QD and CC are equal. The specification 
of this class and the corresponding conditions are completely general and 
apply to any quantum system in a state in this class and satisfying these 
conditions. We further find that the relative contributions of QD and CC 
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can be controlled easily by changing the relative magnitudes of Bi and B2. 
Finally, we connect our results with the monogamy relations between the 
EOF, classical correlations and the quantum discord of two qubits and the 
environment. 

Our work reported in chapter 5 solves the problem of finding a geometric 
measure of QD in a bipartite state. 

As for the future work with quantum discord, we think that finding a 
viable relation between quantum discord and entanglement may be very 
useful for mixed states. Such relations involving inequalities have been ob- 



tained fl36j . Quantifying QD for many particle systems and relating it to 
various physical properties (like order parameters for the case of quantum 
phase transition) is an area we find interesting. We feel that in this area 
QD may be more useful than entanglement because QD measures total 
quantum correlations while in most cases entanglement measures quantum 
correlations breaking Bell inequalities. 
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